Energy distribution of individual quasars from far-UV 
to X-rays: I. Intrinsic UV hardness and dust opacities 



Sinhue A. R. Haro-Corzo^'^, Luc Binette^, Yair Krongold^, Erika Benitez^, 

Andrew Humphrey^, Fabrizio Nicastro^'^ 

and 

Mario Rodriguez-Martmez^ 



ABSTRACT 

Using Chandra and HST archival data, we have studied the individual Spec- 
tral Energy Distribution (sed) of 11 quasars, redshifts 0.3 < z < 1.8. All UV 
spectra show a spectral break around 1100 A. 5 X-ray spectra show the presence 
of a "soft excess" and 7 spectra showed an intrinsic absorption. We found that 
for most quasars a simple extrapolation of the far-UV powerlaw into the X-ray 
domain generally lies below the X-ray data and that the big blue bump and the 
soft X-ray excess do not share a common physical origin. We explore the issue of 
whether the observed SED might be dust absorbed in the far and near-UV. We 
fit the UV break, assuming a powerlaw that is absorbed by cubic nanodiamond 
dust grains. We then explore the possibility of a universal SED (with a unique 
spectral index) by including further absorption from SMC-like extinction. Using 
this approach, satisfactory fits to the spectra can be obtained. The hydrogen 
column densities required by either nanodiamonds or amorphous dust models 
are all consistent, except for one object, with the columns deduced by our X-ray 
analysis, provided that the C depletion is ~ 0.6. Because dust absorption implies 
a flux recovery in the extreme UV (< 700 A), our modeling opens the possibility 
that the intrinsic quasar SED is much harder and more luminous in the extreme 
UV than inferred from the near-UV data, as required by photoionization models 
of the broad emission line region. We conclude that the intrinsic UV SED must 
undergo a sharp turn-over before the X-ray domain. 

Subject headings: quasars: general — X-rays: galaxies — UV: galaxies — galaxies:ISM, 
dust, extinction, nanodiamond — AGN: individua( Pks 1354+19, 3C454.3,Pks 1127-14 
Pks 1136-13, Pks 0405-123, 3C351 PG 1634+706, PG 1115+080, 3C334,B2 0827+24, 
01363) 
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1. Introduction 

The spectra of quasars and Seyfert 
galaxies sliow strong emission lines super- 
imposed onto a bright continuum. The 
continuum contains a significant feature 
in the optical-ultraviolet region, known as 
"the Big Blue Bump" (BBB). The emis- 
sion lines are generally believed to result 
from photoionization, in which the far-UV 
part of the BBB is reprocessed into line 
emission. The shape of the (ionizing) con- 
tinuum in the far-UV domain is therefore 
crucial in the modeling of the observed 
intensities of the emission lines. The con- 
tinuum of Active Galactic Nuclei (AGN) is 
furthermore believed to be responsible for 
the ionization of gas components observed 
in absorption, either within the AGN envi- 
ronment, the so-called "warm absorbers" 



(WA) (e.g. lMathur et al.lll994bl ). or within 
intergalactic space (e.g. Lya forest and Ly- 
man limit systems). Due to the huge pho- 
toelectric opacity of the Galaxy, the Spec- 
tral Energy Distribution (sed) of the ion- 
izing continuum, between the Lyman limit 
and the soft X-rays (EUV-X), is poorly 
known. Fortunately, owing to the redshift 
effect and the transparency of the quasar 
environment, it has been possible to infer 
the SED of quasars down to ~ 300 A. From 
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the X-ray side, the quasar sed becomes 
visible again. The physical relationship 
between the soft X-ray s and the BB B is 
still a matter of debate (jwilkesi [20041 ) . 

Using the same HST-FOS database as 



Telfer et all fl200l TZ02), we present a 
detailed study of eleven quasars for which 
a good quality far-UV spectra exist as well 
as archived Chandra ACIS-S data. Our 
objective in this is to set constraints on 
the SED behavior in the unknown spectral 
domain between the UV and the X-rays. 
Any information about this spectral re- 
gion, even if indirectly gathered, could im- 
prove our understanding of accretion disk 
physics as well as of the physical conditions 
of the broad emission line region (BELR) 
and warm absorbers. 

We now briefly review the salient points 
on what is known about quasar SED, in 
the near and far-UV as well as in the X- 
ray domain. We will explore a novel role 
for dust in the paradigm that the BBB is 
produced by an accretion disk. 

1.1. The far-UV photon crisis 

Satellite and ground-based observations 
of distant quasars showed that the BBB 
peaks in vF,, aroun d 1000 A (rest- f rame) 



flObrien et al.l Il988l : iFrancis et al.l Il991 



Zheng et al. 19971 . hereafter Z97). To char- 
acterize the general trends of the con- 
tinuum shape, various authors co-added 
spectra from samples of quasars at dif- 
ferent redshifts. In Fig.[T] we plot the 
Sloan Digital Sky Survey (SDSS) com- 
posite SED (gray continuous line) con- 
structed by Vanden Berk et aD (2001), us- 
ing 2200 AGN spectra of redshifts between 
0.044 and 4.8. TZ02 showed that the 
sharp turnover shortward of ~ 920 A is 
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mostly due to the cumulative effect of un- 
resolved intergalactic Lya forest absorbers, 
which r esults in the so-called Lym an valley 
trough (IMoller fc .Takobsenl[l99nl ). To cor- 
rect this absorption, TZ02 (and previously 
Z97) applied a statistical correction to each 
spectrum before averaging them. The re- 
sulting composite, represented by the red 
continuous line in Fig.dl consists of 332 
HST-FOS archived spectra of 184 quasars 
with redshifts between 0.33 and 3.6. In 
this case, the resulting SED extends into 
the far-UV, as one would expect if this con- 
tinuum was responsible for the photoion- 
ization of the BELR. A significant break 
occurs around 1100 A, followed by a steep 
flux decrease in the far-UV, equiyalent to 



V 



-1.7 



oc 



Korista et all fll997bD 



pointed out the difficulties in reproducing 
the equivalent widths of the high ioniza- 
tion lines of HellA1640, ClvAA1549 and 
O VI A 103 5, using a power law as soft as 
Note that the TZ02 composite with 
v~^-'^ is somewhat harder than assumed by 
the aforementioned authors. 

A puzzling aspect is the signiflcant vari- 
ation in the hardness of the far-UV con- 
tinuum among different quasars. Many 
quasar spectra reveal an SED significantly 
steeper than the composite, yet there is lit- 
tle evidence that the emission lines are very 
different in those objects. For instance, 
the most extreme case is Ton 34, which 
shows an ultrasteep continuum oc z/~^-^ but 
a near-UV line spectr um not very different 



from other quasars (IBinette fc Krongold 
20071. hereafter BKOT)^ iBinette et al 



(120051, BOS) proposed that crystalline car- 
bon dust, that is, nanodiamond^, might 



^B05 combined two flavors of nanodiamonds, mete- 
oritic (i.e. surface-hydrogenated) and cubic (with- 



be responsible for the 1100 A break, open- 
ing the possibility that a rise in flux may 
take place at energies higher than 20 eV. 
According to this hypothesis, the absorp- 
tion trough would indent the continuum 
only within a relatively narrow spectral re- 
gion, between 10 and 15 eV, followed by a 
flux rise at higher energies. Hence, such 
absorbed SED would still produce more 
photons in the extreme UV beyond 20 eV 
than obtained by extrapolating the TZ02 
composite. Examples of SEDs absorbed by 
nanodiamond dust are represented by the 
two dotted lines in Fig.[T|, corresponding to 
Models I and II, both deflned in §[6l The 
sugge stion that AGN m ay be a "factory" of 
dust (lElvis et al.ll2002l) or nan odiamonds is 
not new (iRouan et aDl2004al Jbll 



Korista et al.l (Il997bl ) speculated that 



there could be a second bump peaking 
near 54 eV that would provide the nec- 
essary hard photon luminosity. In any 
event, state-of-the art photoionization cal- 
culations of the BELR show a preference 
for a much harder ionizing SED compared 
to what is commonly observed in quasars. 
To illustrate this, in Fig. [J we show t wo dis- 



tributions used by iKorista et al.l (Il997al . 



K097) in their extensive grid of BELR 
models: the short and long-dashed cyan 
lines correspond to an exponential cut-off 
SED of the type u'^'^ exp—hu/kT, with 
kT = 43 and 86 eV, respectively. The Lo- 
cally Optimall y emit ting Cloud model of 



of Casebeer et al. 



Baldwin et al.l (11995 ) and, more recently. 



use a similar 



SEe|^. Also shown in the same figure is the 



out surface adsorbates), assuming the small size 
regime. The resulting extinction curves Al and 
Dl are reproduced in Fig. [HI 

^Although Casebeer et al.l ( 20061 ) find evidence of 
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popul ar quasar SED of iMathews fc Ferland 
( 119871 ) (yellow dashed line). This inconsis- 
tency between the SED preferred in pho- 
toionization calculations and the SED ac- 
tually observed in most quasars is unfor- 
tunate and requires a solution. 

An additional aspect to be taken into 
account is that nearby, less luminous AGN 
might possess a n intrinsically harder SED, 
as reported by IScott et al.l (120041 ). who 
used FUSE data to construct a composite 
SED comprising AGN of redshifts < 0.7. 
These authors did not find evidence of a 
UV-break, at least for an important frac- 
tion of their objects. The 11 quasars re- 
ported in this paper, however, are consis- 
tent with the presence of a break. As in 
BOS, our analysis will focus on individual 
spectra rather than on composites. 

1.2. Constraints provided by the X- 
ray domain 

Since the X-ray domain undergoes lim- 
ited absorption as compared to the Ly- 
man continuum, it provides us with con- 
straints on how the BBB might behave 
in the unobservable extreme-UV domain. 
This unknown domain is depicted by the 
shaded area in Fig.[Tl Past studies (e.g. 
Reeves &: Turnerl l2000l . and references 
therein) revealed that the typical spec- 
trum of a quasar in the hard X-ray band 
(> 2keV) is dominated by a powerlaw- 
like emission flux of index ax ~ — 1- It 
has also been found that a gradual up- 
turn occurs below 2keV, the so-called X- 



ray soft excess (jComastri et al.lll992l ). The 



a very hard SED in RE 1034+39, they also report 
that "the slope appears to become flatter in the 
FUSE spectrum" shortward of Ovi, which they 
qualify as " an unexplained feature" that is left out 
of their analysis. 



latest generation of X-ray satellites (XMM- 
Newton and Chandra), owing to their im- 
proved spectral resolution and sensitivity. 



Gambill et al. 


2003: 


Porquet et al. 


Piconcelli et al. 


2005 


, hereafter 



However, its nature as well as its relation- 
ship with the BBB are still a matter of 
debate. 

1.2.1. The broken powerlaw X-ray model 
Considerin g a sample of 23 lo w redshift 



quasars [the iLaor et al.l (119971 ) subsam- 



ple of P G quasars] obsery e d wit h XMM 



Newton, iBrocksopp et al.l (120061 ) favored 



a broken powerlaw model to fit both the 
soft energy excess and the harder X-ray 
segment. They found that the two indices 
used in their models are well correlated, 
and they concluded that the X-ray com- 
ponents in the range 0.3-10 keV shared a 
common origin. These authors also re- 
ported that accretion disk parameters re- 
sponsible for the UV BBB appeared to be 
independent of the 0.3 keV excess. 

We averaged the brok en powerlaw fits 



to 13 AGN {z < 0.4) that IBrocksopp et al. 



( 120061 ) reported as 'good fits' (see their Ta- 
ble 8). Such average is represented by the 
two navy blue lines in FigJTl Each line illus- 
trates a predefined valuci of the aox index 
(-1.3 and -1.5). The flux from the TZ02 
composit^ at ~ 2500 A was used to de- 



^Negative indices imply decreasing flux from the 
UV to the X-rays. 

''Using the observed composite flux rather than 
that of a continuum model is preferable since it 
includes the contribution from the 'small bump' 
around 2700 A, which is believe d to be due to un- 
resolved multiplet Fen emission ( Wills et al.lll985[ 
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termine this index. iLaor et al.l (119971 ) and 
TZ02 suggested that the soft excess might 
correspond to an extension of the BBB 
far-UV powerlaw into the X-rays. Our 
Model II in Fig.[Tl for instance, would con- 
nect smoothly with the soft excess if ex- 
trapolated beyond the shaded area. Al- 
though a tempting proposition, we should 
beware of general conclusions derived from 
composite SEDs. A detailed analysis cen- 
tered on individual quasar SEDs leads us to 
reject the above premiss in §|H 

The SED behavior within the data gap 
(the shaded area in Fig.[l]) is not known 
a priori. In order to establish useful con- 
straints on the different BBB models, we 
will extrapolate the UV models across this 
gap and compare them with the X-ray ob- 
servations. 

1.2.2. A thermal model for the X-ray ex- 
cess 

Using a larger subset of 40 PG (Palomar- 
Green) quasars, PI05 fit five different mod- 
els to the soft X-rays: either one or two 
blackbodies (A, E), a multicolor blackbody 
(B), a bremsstrahlung emission model (C), 
and a powerlaw (D). Except for one object, 
all these models require the addition of an 
underlying powerlaw in order to fit the 
high energy band. Models A or E were 
favored for 19 objects, while model D was 
favored for 6 objects only. The soft ex- 
cess in some cases shows a top-hat behav- 
ior, which favors a thermal interpretation. 
The green line in Fig. [1] represents an aver- 
age of 13 model fits (blackbody+powerlaw) 
of quasars with z < 0.4 and of type A 
and E, using parameters of Table 5 from 



PI05. This average SED has been arbitrar- 
ily scaled to an aox of —1.45, and the mean 
blackbody temperature describing the ex- 
cess is ~ 140 eV. PI05 emphasize that the 
five functional forms they used must be 
considered as phenomenological and not 
as physically consistent models. 

A well resolved rollover (i.e. decreas- 
ing fiux towards the UV) of the soft X- 
ray component, or at least a strong cur- 
vature, has also been reported in Ark 120 
fIVaughan e t al.ll2004h Mrk359 dO'Brien et al, 
200l|) and NGC4052 dCoUinge et al.ll200lh . 
In a similar fashion to PI05, we adopt in 
this Paper a single blackbody to fit the soft 
X-ray excess. 

1.3. The physical origin of the BBB 

Traditionally, the BBB has been in- 



disk 


emission ( 


Shields et al. 


1998: 


Malkan 


1983 




Edelson & Malkan 


19861) around a 



supermassive blackhole (BH). The pres- 
ence of a BH in the nucleus of all grand de- 
sign galaxies has since bee n confirmed (e.g. 



Ferrarese &: Merritt 2000, and references 



therein). A persisting difficulty, however, 
is that in the region of the far-UV break, 
the BBB is not fitted well by 'bare' accre- 
tion disk models. This is illustrated by a 
comparison of the composite SED from Z97 
with sta te-of-the-art di s k mo dels such as 
those of iHubeny et al. fl2000h (c .f. their 



Fig. 22), which include H and He opac- 
ities and non-LTE transfer. A solution 



Yip et al.ll2004f) . 



proposed by Z97 and IShang et al.l (12005 
SBG05) is the comptonization of the disk 
emission by a hot corona above the disk. 
This corona would smear the intrinsic disk 
features such as the Lyman break and re- 
sult in a SED with a powerlaw tail in the 
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far-UV. Hence accretion disks do not nec- 
essarily imply a sharp thermal cut-off, as 
predicted by bare disk emission models. 

For a given SED optical flux, the effi- 
ciency in generating ionizing photons is 
higher for hotter accretion disks, and these 
imply a fairly hard continuum in the near- 
UV, wit h indices aAf riv > -0.2. The TZ02 



and the IScott et al.l (120041 ) composites, on 
the other hand, reveal a much softer con- 
tinuum with a mean a^uv index of —0.67 
and —0.57, respectively. In this Paper, we 
explore whether dust might be at the origin 
of this apparent SED softne ss in the near- 



UV. F ollowing the results of iHopkins et al. 
( 12004J ). who finds evidence of extinction in 
AGN similar to that encountered in the 
Small Magellanic Cloud (SMC), we will be 
considering such extinction in §[6] to dered- 
den the observed SEDs until they match a 
target spectral index. 

As for the X-rays, a popular physical 
scenario for the hardX-raj (> 2 keV) com- 
ponent is that of a hot corona placed above 
the accretion disk, which comptonizes the 
UV to soft X-ray thermal emission from 
the disk a nd up-scatters it into the ha rd X- 
ray band (IHaardt fc Maraschil Il99ll ) . On 
the other hand there are many competing 
physical mechanisms that have been pro- 
posed to explain t he soft excess. They are 
revi ewed in PI05, iBrocksopp et al.l (120061 ) 
and IPorquet et al.l (1200 



1.4. Objectives of this work 

The issues that will be tackled in this 
Paper are the following: using individual 
spectra rather than composites, we ana- 
lyze to what extent the extrapolation of the 
far-UV results in an X-ray flux compatible 
with Chandra observations. In the event 



that the UV break was the result of crys- 
talline dust absorption, we wish to know 
whether the rollover proposed by BOS at 
18.5 eV, followed by a powerlaw tail in the 
extreme-UV, is compatible with the ob- 
served X-ray flux. Is this powerlaw tail an 
acceptable function? Or is a thermal cut- 
off more natural? Is SMC-type dust able 
to account for the observed wide range in 
near-UV spectral indices? Are the columns 
of material implied by all the dust com- 
ponents compatible with the intrinsic ab- 
sorption column inferred from the X-rays? 
What is the impact of the dust-absorbed 
SEDs studied in this Paper on the emis- 
sion spectrum of the Narrow Line Region 
(NLR)? The framework that we introduce 
to explore these issues should be useful for 
follow up studies concerned with the far- 
UV spectral gap in AGN SEDs. 

2. The observational dataset 

In this Paper, we study 11 quasars for 
which high quality spectral data sets ex- 
ist in both the HST-FOS and the Chan- 
dra archives. Of these, 9 are radio-loud 
quasars (RLQ) and 2 radio-quiet quasars 
(RQQ). The larger fraction of RLQs is 
due to small number statistics of quasars 
common to both databases. In Tabledl 
we list the object names (Col. 1) in or- 
der of decreasing X-ray counts. The italic 
letter shown in Col. 2 is used for identi- 
fication throughout the Paper and in all 
figures. Other table entries are the red- 
shift (Col. 3), the Galactic Hi column den- 
sity (iDickev fc Lockmanlll99oh (Col. 4), ra- 
dio loudness of each source, RQQ vs RLQ 
(Col. 5), and the UV spectral class accord- 
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ing to the nomenclatur^ of BOS (Col. 6). 
The X-ray journal of observations is also 
briefly described and includes the follow- 
ing: the identification number for each 
source in the Chandra database (Col. 7), 
the exposure time (ksec) of each observa- 
tion (Col. 8), the total number of counts for 
each spectrum (Col. 9), the frame-time in- 
formation (Col. 10), and Chandra's obser- 
vation date (Col. 11). HST's observation 
date is also listed (Col. 12). 

We have postponed the analysis of the 
14 quasars (10 RLQ and 4 RQQ) that 
showed evidence of a WA in the UV to a 
future paper, except for 3C351, for which 
the detailed analysis is presented in Ap- 
pendix|Xl Since it is the only object in 
common with the SBG05 sample, its inclu- 
sion is intended to facilitate a comparison 
of our results with theirs. 

2.1. The UV Data 

The UV spectra for our sample were 
extracted from the TZ02 database, which 
R. C. Telfer was kind enough to provide. 
It comprises spectra, mostly HST-FOd§, 
of 184 quasars, already reduced, corrected 
for Galactic dust extinction, de-redshifted 
and finally rebinned uniformly. Further- 
more, the spectra have been corrected by 
TZ02 for the presence of Lyman limit ab- 
sorbers (down to r > 0.3), as well as 
for the Lya absorption valley (caused by 
the cumulative absorption from unresolved 



^ Class A spectra show a UV break similar to that 
seen in the composite SED of TZ02, while class B 
quasars show a much steeper flux drop shortward 
of the break (see classification scheme in Fig 2 of 
BOS). 

^Note that the TZ02 sample includes 3 HST-STIS 
and 6 HST-GHRS spectra. 



Lya forest lines, see § 11. ip . We selected 
spectra that possess the following charac- 
teristics: an underlying continuum of ac- 
ceptable S/N (see Col. 3 in TableEl suf- 
ficient spectral coverage above and below 
the 1100 A UV break (see TableE]), and the 
absence of deep absorption lines or similar 
features from a BALQSO. 

The rest-frame UV spectra and their er- 
ror bars are plotted in loguF^, vs. logu 
in FigS To increase the S/N and to re- 
duce the smearing from undifferentiated 
high amplitude scatter, these spectra have 
been rebinned to group together 10 points 
per resolution element. 

2.2. The X-ray Data. 

Our sample includes quasars for which 
Chandra's observations were carried out 
with ACIS-S. The corresponding datasets 
were freely available from the Chandra X- 
Ray Centeilll. 

The files were processed using the CIAO 
software (version 3.0.2) and following the 
on-lin^data analysis 'threads' provided by 
the Chandra X-ray Center in order to ex- 
tract source and background spectra. Re- 
distribution Matrix Files (RMF) and Aux- 
iliary Response Files (ARF) were gener- 
ated by using the CIAO tools mkrmf and 
mkarf, respectively. 

We estimated whether each source was 
affected by pile-up (at a level > 10%) us- 
ing the expression PU = where 
TC is the Total Counts, FT is the frame- 
time and ET is the Exposure-Time of each 
observation (see Tabled]). If PU > 0.3, 
a pile- up correction was carried out. We 



''http://chandra.harvard.edu/ 

^http: / /cxc. harvard.edu/ciao /threads / data.html 
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found that only quasars b, c, g and h 
showed a significant pile- up. For these, 
the correction was performed by select- 
ing only counts within an annulus cen- 
tered on the source and defined by radii 
0.6-3 arcsec. We took into account that 
a proper fit requires at least 1000 counts. 
For all other objects, we used a circle of 
3 arcsec radius centered on the source. The 
selection of this size aperture guarantees 
that counts due to extended radio-jets, if 
present, are avoideci^ as well as counts 
from weaker nearby sources. This pro- 
cedure was performed by using either a 
box or a circle around each non-nuclear X- 
ray source. For each source, we generated 
a background spectrum by considering all 
the counts within an annulus defined by 
the radii of 10 and 20 arcsec, again cen- 
tered on the object. 

3. Separate model-fitting of the UV 
and X-ray segments 

We use powerlaw segments to char- 
acterize the behavior of the energy dis- 
tribution in the UV and in the X-rays. 
We performed spectral fitting^ of the 
11 quasars, using the Sherpallj package 
in CIAO. Each powerlaw segment is ex- 
pressed by the form Fj 
space Fi 



(erg cm ^ s ^ Hz 



1 



'^'^ in flux 
or by the 



®We had no means to remove the contribution 
from the unresolved X-ray jet component close to 
the nucleus, whic h is possibly amplified due to 
Doppler-beaming (jGambill et a l. 2003). The X- 
rays measured in RLQs might therefore be con- 
sidered upper limits of the nuclear fluxes. 
^•^The commands used for the 
different fits are described in 



form uFi, oc in the customary uF^, 

representation. Since most figures are in 
log uF^ vs. log u, we prefer to use the in- 
dex throughout most of the text. Such an 
index directly reflects the apparent slopes 
of the powerlaw segments seen in the fig- 
ures. We recall that ai, = [3^ — 1. A few 
authors prefer to use an index defined in 
Fx space {Fx oc A"*"^^), in which case, we 
have 7a = — 2 — a,^ = — 1 — For the X- 
rays, the photon index Fx commonly used 

in the literature is given by Fx = 2 — /5x 
F V 

{-j^ oc ^, note sign convention for Fx). 

3.1. Modeling of the ultraviolet re- 
gion 

As in TZ02, we fitted two broken pow- 
erlaws, one for the near-UV and the other 
for the far-UV, of indices /3^uv ^'^'^ /^fiiv 



respectiveljo. These preliminary fits as- 
sume that no intrinsic dust is present. We 
employ the superscript 'dfe ' to label quan- 
tities related to the 'dust-free' fits. The 
wavelength at which the two broken pow- 
erlaws join is defined as A^^^. 

Spectral fitting was carried out in two 
different ways: the first technique con- 
sisted of fitting a broken powerlaw (bp) to 
all the line-free continuum regions, leav- 
ing out the important emissi on lines, in 



a similar fashion as TZ02 and Scott et al. 



(I2OO4J ). This masking is equivalent to ex- 
cluding the following wavelength intervals 
from the fit: 750-800, 820-850, 900-1095, 
1100-1135, 1150-1265, 1380-1420, 1470- 
1610, 1830-1950, 2700-2880 A0. The sec- 



http : //cxc . harvard . edu/ ciao/ ahelp/xs . html 

""^""^http:/ /cxc. harvard.edu/ciao/ahelp/sherpa. html 



^^Note that TZ02 used indices, and that = 
-f 1 

^^Objects e, g, h, i and k required slight modifica- 
tions of these blackout intervals, c.f. App.El 
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ond technique consisted of fitting the con- 
tinuum and the emission hues together, 
using a hne profile consisting of both a 
narrow and a broad gaussian component. 
We found good agreement between the 
spectral indices derived using each of the 
two techniques, and therefore we consider 
both reliable. Since we are primarily in- 
terested in the continuum rather than in 
the lines, we adopted the masking tech- 
nique due to its simplicity. Because the 
continuum break occurs near 1100 A, that 
is, within the 1100-1135 mask or the Lya 
mask, the determination of A^^^ was com- 
plicated and required an initial eye-ball es- 
timate. The position of A^^^ was deter- 
mined from an initial fit that did not con- 
sider any masking. We then froze Xfjl and 
proceeded to determine new values of /3^uv 
and Pp^y by least-square-fitting within the 
unmasked regions only. These values are 
presented in TableO The two segments 
from the best fitted broken powerlaw are 
superimposed onto the UV spectra in Fig. [2] 
a-k. The far-UV break near llOOA IS ap- 
parent in each quasar spectrum. The mean 
break position characterizing our sample is 
Xf^l = 1165 A, which is significantly red- 
ward of the Lyman limit. 

The index difference, /3f{jv~/^nijv; gives 
us a measure of the steepening across the 
break. The values we determined range 
from —1.93 to —0.72, with a mean differ- 
ence value of —1.39. The latter is compara- 
ble to the mean index difference of —1.22 
and -0.81 for the RLQ and RQQ of the 
TZ02 sample, respectively. 

3.2. Modeling the X-ray region 

To model the intrinsic X-ray continuum 
of each quasar, the data were first rebinned 



in such a way that the number of counts in 
each bin exceeded 20. This had the ad- 
vantage of preserving the Poisson statis- 
tics when fitting Gaussian profiles and of 
facilitating the statistical handling of the 
background substraction. The observed 
continua with their error bars are shown 
as counts s~^ keV~^ vs. energy (observer- 
frame) in the top panels of Figs. [3] a-k. 

The X-ray spectra were fitted in the 
observer-frame by models that were folded 
into the spectral characteristics of the 
ACIS-S instrument. We initially fitted a 
simple powerlaw attenuated by Galactic 
absorption of column N§^. This initial fit 
was limited to the range 2.5-6 keV (rest- 
frame), thereby excluding the soft X-ray 
region and the FeKa line region around 
6.4 keV. It did not consider possible ab- 
sorption by moderate columns of neutral 
gas (i.e. < 10^^ cm~^) or by larger columns 
of ionized gas (i.e. < l O^^cm"^), both 



intrin sic to the quasar (IKrongold et al. 



20031 ). The residuals of this initial fit 



are presented in the middle subpanels of 
Figs. [3] a-k. 

For almost all the quasars, such a simple 
model could not provide a satisfactory fit 
(see Table[3]) over the entire energy range 
(i.e. 0.3-6 keV observer- frame) . The mid- 
dle subpanels show positive as well as neg- 
ative deviations (with a significance > 2cr). 
Shortward of 2.5 keV, we find either a soft 
excess or evidence of absorption gas. Long- 
ward of 4 keV, there is some evidence of a 
line emission feature, which might corre- 
spond to a fiuorescent transition of neu- 
tral iron emission from FeKa at 6.4 keV. 
The next step was to repeat the fit, leaving 
the spectral index of the powerlaw frozen 
to the previous value, but now including 
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a blackbody and/or an intrinsic absorber 
component in order to account for the ob- 
served residuals (see middle panels). We 
left the normalizations of the powerlaw, 
blackbody and intrinsic absorption free to 
vary. We can see in the bottom panels 
a noticeable reduction of the residuals at 
low energies. Note that the fit to 3C351 
(object c) requires a warm absorber com- 
ponent, as shown in App.|A]-c. In this last 
model and for the 2 quasars a and e that do 
not show a clear evidence of soft excess, we 
used an f-test to determine the level of sig- 
nificance between models with and without 
a blackbody. The level of confidence is less 
than 10~^, which confirms that including 
a blackbody component for both quasars 
provides a better fit. 

The continuous line in the top subpanels 
of Fig. [3] represents the final best-fit multi- 
component model, while the correspond- 
ing residuals are shown in the bottom sub- 
panels of Fig. [31 In the particular case of 
PG 1115-1-080, we had to consider a broken 
powerlaw that joins at 1.69 keV (see note 
k in App.Rl). TableO lists all the param- 
eters of the final model. Uncertainties at 
a 2(T-confidence level regarding the power- 
law index of all quasars were derived, con- 
sidering only the 2.5-6 keV (rest-frame) 
range. 

To summarize, the inferred X-ray spec- 
tral energy distribution of each quasar in- 
cludes the following components: thermal 
emission from a blackbody (zbb) for ob- 
jects a, e, f, g, i, and a Gaussian Profile (g) 
representing the FeKa line for objects c, g, 
i, j. Evidence was found for neutral gas ab- 
sorption (of column N^) intrinsic to the 6 
objects a, b, d, e, f, h and a WA in the ob- 
ject c, while for the remaining 4 objects, we 



froze the other parameters at their previ- 
ous best-fit values and calculated an upper 
limit for based on a confidence level of 
2a. 

4. Combining the UV and X-ray 
SEDs in the dust-free case 

In this Section, we assume the environ- 
ment of quasar to be dust-free and make 
use of a broken powerlaw to describe the 
UV spectrum of each quasar. 

4.1. Matching the UV and X-ray 
components 

We used IDL v. 5.5 to combine the UV 
and X-ray segments. In FigjH we over- 
lay in uFiy the UV and X-ray model fits 
to obtain a single tentative description of 
each SED, as a function of u (rest-frame). 
In order to avoid any clutter due to undif- 
ferentiated scatter and overcrowding, the 
original data from Figs.[3] and H] have been 
left out. Arrows near the top of each panel 
delineates the spectral extent of the origi- 
nal data. 

The near- and far-UV powerlaw fits are 
shown as dashed lines, while the two cross- 
ing dotted lines represent a 2a deviation 
about the best-fits. We recall that the in- 
dices Pi, (with a,^+ 1) are explicitly 
listed in each panel and represent the ac- 
tual sloped in uF^, plots. 

As for the X-ray domain, the previous 
best-fit models are plotted as dashed lines. 
The associated limit variations, assuming 
a 2a uncertainty for the parameters used, 
are represented by dot-dashed lines. These 



The conversion from other index definitions such 
as Fx and 7a is obtained as follows: /3x = 2 — Fx 

and Pi, = -7a - 1- 
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'fluxed' X-ray models represent the intrin- 
sic SED corrected, not only for Galactic ab- 
sorption (GA) by column N§^, but also 
for intrinsic absorption by gas column N§ 
(§ l3.2p . A vertical dotted- line denotes the 
lower energy boundary of the data. 

We determined the energy, Ej-^J^, of the 
intercept that occurs when extrapolating 
both the far-UV and X-ray segments. In 
Fig.m a small circle denotes the position 
of the intercept. The mean intercept en- 
ergy for the sample is Ej{^^^ = 265 eV (this 
average excludes two objects without a de- 
fined intercept, as discussed below). The 
symbol '@' represents the UV break po- 
sition Xf^^f^. We also derived the aox in- 
dex, that is, the spectral index of an imag- 
inary powerlaw joining the two flux val- 
ues at 2500 A and 2keV (as denoted by 
open triangles). This index is given in Ta- 
ble. H] and in the inset of Fig.H] (along with 
\fj.^f^ and Ejl'^^). Its value ranges between 
— 1.78 and —1.11, with a mean of —1.46. 
The aox values from samples of QSOs of 
comparable luminosit ies tend to fall in the 
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4.2. Discussion I: the dust-free case 

We now discuss the issue of whether the 
extrapolation of the far-UV powerlaw leads 
to a flux level compatible with the one in- 
ferred from the X-ray model. There are 
two relevant cases to consider: A- the ex- 
trapolated far-UV flux model lies below or 
at the level inferred from the extrapola- 
tion of the X-ray model, B- the extrap- 
olated far-UV lies above the level inferred 
from the extrapolation of the fluxed X-ray 
model. In case A, at a phenomenological 



level, the extrapolation leads to a congru- 
ent result, that is, the extreme-UV flux 
either vanishes relative to the soft X-ray 
component or connects smoothly with it, 
while in caseB, we conclude that a sim- 
ple extrapolation of the UV powerlaw is 
nonsensical. The far-UV segment in case B 
objects must instead undergo an abrupt 
steepening before reaching the X-ray do- 
main. With case A, the extrapolation of 
the far-UV segment does not lead to in- 
consistencies with the X-ray data, while it 
does with caseB, unless further steepening 
occurs in the extreme-UV, possibly in the 
form of an exponential turn-over. 

The statistics are the following: out of 
11 quasars, 9 are case A and only 2 caseB 
[b, k). We conclude that the extrapola- 
tion of the far-UV powerlaw leads to a flux 
level compatible with that measured in the 
X-rays, except for two sources. This does 
not entail that the continuum in the unob- 
served region between the far-UV and the 
X-rays necessarily behaves as the partic- 
ular functional form considered. The in- 
trinsic SED may still steepen at some un- 
known intermediate energy (between 40 
and 350 eV), before the soft X-ray excess 
component takes over. Such a possibility is 
favored in scenarios in which a thermal dis- 
tribution is believed to provide a more ap- 
propriate description of the BBB. Finally, 
absorption by dust can strongly affect the 
observed SED in the far-UV. Correcting for 
such absorption would lead to a different 
intrinsic SED and its extrapolation may re- 
sult in a different balance between case A 
and caseB, as is found in §|6l 

Quasar variability is a concern, since the 
observations in the UV and X-rays were 
not simultaneous. The characteristic am- 
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with source luminosity ( 


Pica & Smith 


1983: 


Trevese et al. 


1994; 
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ook et al.lll994 


: iHawkins 


1996 


: Cristiani et al. 


1997). It is smaller 



for qu asars than Seyferts. iGiveon et al. 



(119991 ) found a B band rms variability 
of 0.14 mag, following a 7 year monitor- 
ing campaign of 42 PG quasars of redshift 
< 0.4. Eyen though such amplitude is rel- 
atively small, there are clear indications 
that the amplitude o f variations in the far - 
UV are larger (e.g. iReimers et al.l l2005l ). 
As for the X-ray domain, it has been shown 
in general terms that va riability anti- 



corre lates with luminosity (INandra et al. 



19971 ). and our sample consists of high lu- 
minosit y quasars. For our sa mple in par- 
ticular, iGambill et al.l (120031) reported no 
short-term high-amplitude variability for 
the quasars b, /, and i. As for longer term 
variability, these authors reported a de- 
crease of a factor ~ 1.5 over nine years 
for quasar b and a factor ~ 3 decrease over 
twenty years for quasar /. It is conceivable 
that when the X-ray and UV observations 
are not concomitant, the distinction be- 
tween case A and caseB for a particular 
quasar may depend on the temporal lim- 
itations of the data. On the other hand, 
when a reasonably sized sample is consid- 
ered, flux variations should have a neutral 
impact on the relative proportion of case 
A vs. B. We note that the ratio of the ex- 
trapolated far-UV model to the observed 
X-ray flux, both estimated at the long 
wavelength limit of the X-ray data, exceeds 
±0.6 dex for 7 of the 11 quasars. We there- 
fore consider it unlikely that our results of 
a majority of caseA's could be reverted 
owing to AGN variability. The possibility 
that the X-ray continuum is significantly 
contaminated by an unresolved jet com- 



ponent, however, might have caused real 
case B's to be misclassified as case A's. 
This potential obstacle cannot be removed 
without higher spatial resolution and un- 
til the issue of the real nat ure of the soft 
excess becomes clearer (see Gambill et al. 



20031 . and references therein). 

As an exercise, we compared our re- 
sults with those of SBG05, who built com- 
prehensive SEDs of nearby AGN that ex- 
tend from the IR up to the X-rays. Their 
X-ray data were taken from the ROSAT 
'AH Sky Survey' and public PSPC obser- 
vations. Similarly to us, SBG05 fitted a 
powerlaw to the far-UV segment. They 
represented the ROSAT data by a simple 
powerlaw. Among the 16 AGN, for which 
they had X-ray data, we determined from 
their Fig. 3 that 6 are case A and 10 are 
caseB. We do not have a simple explana- 
tion for the higher frequency of case B's 
in the SBG05 sample. One possible ex- 
planation is that the soft X-ray excess is 
not satisfactorily accounted for by a pow- 
erlaw alone. We recall that we added a 
blackbody to the underlying powerlaw for 
our sample, in order to fit the excess, and 
that we found a noticeable soft X-ray bump 
in 5 objects (§[3]). Alternatively, since the 
SBG05 sample consists of nearby (less lu- 
minous) AGN, systematic differences in 
the energy distr ibution could exist (see 



Scott et al.l I2OOJ) that would result in a 
larger proportion of case B's in SBG05. In- 
terestingly, we have one object in common 
with the SBG05 sample, 3C351 (labeled c), 
which according to their Fig. 3 is caseB, 
while it turns out to be case A in our anal- 
ysis. A thermal bump appears not to be 
necessary to model this quasar, but we find 
evidence of a strong WA, which obviously 
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somewhat complicates a direct comparison 
with SBG05. The WA model included in 
our fit is described in AppendixRl 

If we consider separately the RQQ or 
the RLQ composite SEDs from TZ02, they 
both fall in the case A category. There- 
fore, regarding the dominance of case As 
in the dust-free case, composite SEDs as 
well as individually analyzed spectra lead 
to similar conclusions. Case A does not en- 
tail, however, that the X-ray and far-UV 
components share the same physical origin. 
This was conjectured in previous works 



fe.g. iPuchnarewicz et al.lll996l : lLaor et al. 



19971 ). in which the soft excess was con- 
sidered to be the prolongation of the far- 
UV component observed in the compos- 
ite SED. The picture that is emerging 
from recent works rather suggests the op- 
posit e. For instance, using X MM-Newton 
data, iBrocksopp et al.l (120061 ) and PI05 ar- 
gue that the soft X-rays are a distinct 
component from the BBB (see discussion 
in § 11.21) . Our analysis supports this in- 
terpretation, since a large fraction of ob- 
jects (5/11) presents a differentiated ther- 
mal bump in the soft X-rays, and none is 
found where the soft X-rays behave as the 
prolongation of the far-UV powerlaw. 

5. Combining the UV and X-ray 
segments in the 'break-corrected' 
case 

5.1. Competing explanations for 
the UV break 

The composite SED inferred by TZ02 
appears to be too soft to account for 
the broad emission lines (see § 11.11 and 
Korista et aTlll997bl l We have initiated 
an ongoing project to explore alternative 



solutions to this problem, assuming that 
the far-UV break is more akin to a lo- 
calized continuum trough, followed by a 
marked recovery in the extreme UV, which 
is the energy region responsible for the high 
excitation emission lines. Various mecha- 
nisms that could ge nerate such a trough 
are summarized by iBinette et al.l (120071 ). 
State-of-the-art calculations of the SED 
for standard geometrically-thin optically- 
thick accretion disks do not reproduce 
the observed break satisfactorily (§ 11. 3p . 
but such models assume a stationary disk 
with a vertical structure supported only 
by gas and radiation pressure. To our 
knowledge, no detailed SED calculations 
assuming a non-stationary disk or with 
an accelerating wind have been carried 
out. Absorption by intergalactic Hi or 
intergalactic dust h ave been discarded by 



Binette et al.l (20031) a nd B05, respectively. 



Eastman et al.l (119831 ) proposed that Hi 



absorption by local clouds accelerated up 
to 0.8c could generate a steepening of the 
transmitted continuum. In order to re- 
produce modern data, this model would 
need to be fine-tuned and extended to the 
extreme UV where the flux recovery is ex- 
pected. More recently, BOS have proposed 
that the break could be the result of ab- 
sorption by crystalhne dust grains local to 
the quasars. 

More work is needed to falsify some 
of these competing explanations. How- 
ever, whatever the physical origin of the 
break, once we can reasonably reproduce 
its shape, we can explore the possibility of 
a universal quasar SED, in which the large 
variations in observed spectral indices are 
the result of absorption by standard dust 
models either SMC-like or Galactic. This 
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possibility of a universal SED is pursued in 
§[11 assuming that the break is due to crys- 
talline carbon dust, as explained in detail 
in the current section. We will assume that 
a SED turnover takes place, not at the ob- 
served break, but at shorter wavelengths, 
around 670 A, as in BOS. 

5.2. Crystalline dust to model the 
far-UV break 

BOS showed that absorption by a novel 
dust component made of crystalline car- 
bon (i.e. nanodiamonds) could reproduce 
the position and detailed shape of the 
UV break in SO objects, out of a total 
sample of 61 quasars with multi-grating 
spectra extending down to at least 900 A 
(rest-frame). An example of such spec- 
trum is that of PG 1148+S49 {z = 0.969) 
shown in Fig.[5l An important predic- 
tion of BOS is that a flux rise should take 
place near 700 A and this rise was ob- 
served in HS 1307+4617, PG 1008+1319 
and Pks 0232-04. BK07 recently pre- 
sented evidence of a far-UV rise in Ton 34, 
the quasar with the most extreme break 
known. The gray dashed-line in Fig. [5] il- 
lustrates the far-UV rise predicted by the 
absorption models of BOS. We now re- 
view specific features of the crystalline dust 
model. 

Turnover in the EUV-X region. If 

the UV-break is the manifestation of dust 
absorption, the intrinsic continuum (once 
dereddened) cannot rise indefinitely in 
pF^, as implied by a single index powerlaw. 
The intrinsic SED must present a turnover 
before the X-ray domain. The aox index 
provides us with loose constraints on such 
a turnover. Typically, aox lies in the range 



—1.4 to —1.6. The radio-loud quasars are 
X-ray stronger, with an aox reaching —1.2. 
This index gives us some clues about the 
drop in flux between the far-UV and the 
soft X-rays, but not on where exactly the 
actual turnover takes place. 

BOS reported evidence of a rollover 
of the continuum in the far-UV SED of 
the well studied high redshift quasars 
HS 1700+6416 and HE 2347-4342 (both 
at Zq ~ 2.8). They parameterized the 
shallow turnover observed by using the 
following multiplicative function: Crou = 
(1 + [Xroii/ ^]~'^^)~^^^ , which produces a 
steepening centered on Xroii- ^ is the pow- 
erlaw index increase and / a form factor. 
The cut-off can be progressive or sharp, 
depending on /. When the SED is multi- 
plied by this function, an index change (6) 
takes place at wavelength Xroii- The values 
inferred from the above two quasars were 
/ = 2.8, Xroii = 670 A (18.SeV) and an in- 
duced steepening of magnitude 6 = —1.6. 
By incorporating this rollover into their 
intrinsic powerlaw SEDs, BOS found a sig- 
nificant improvement in the synthetic sim- 
ulation of the TZ02 composite. Further- 
more, the fit to the far-UV rise observed 
in four quasars was also improved (BOS, 
BK07). We will test whether this func- 
tional form, which behaves as a powerlaw 
in the unobserved domain, is at least con- 
sistent with the X-ray observations. The 
rationale behind this form, rather than a 
thermal/exponential rollover, is the same 
as that proposed by Z97 to model the com- 
posite SED, that is, comptonization of the 
disk emission by a hot corona. 

A simplified nanodiamond dust model 

BOS could reproduce the UV break of 
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quasars, using two kinds of crystalline 
carbon grains: terrestrial cubic and me- 
teoritic. We hereafter adopt a simplified 
version of that model, based on a single 
albeit modified extinction curve consisting 
of terrestrial cubic diamonds only. This 
new curve D3, shown in Fig[6l differs from 
the Dl curve of BOS in that the grain size 
distribution (oc a~^'^) covers a wider size 
range of 3-200 A (instead of 3-25 A). The 
main advantage is that a single extinc- 
tion curve now suffices to model the break 
without the need of meteoritic grain^. 
This is shown in Fig.O where the extinc- 
tion by curve D3 provides a satisfactory 
fit (continuous line) to the sharp break ob- 
served in PG 1148+549. In the same fig- 
ure, the dotted line represents the intrinsic 
SED adopted for this quasar, which con- 
sists of a powerlaw of index /?nuv = 0.8, 
multiplied by the function Crow defined 
above, which produces a shallow rollover 
near 670 A. The new fit assuming model D3 
is quite similar to the one obtained by B05 
(dashed gray line), who used a combination 
of extinction curves Al and Dl. All our 
dust models were computed assuming solar 
metallicity for carbon (C/H= 3.6 x 10~^) 
and full depletion of carbon onto dust (i.e. 
Sc= 1.0). This choice is purely for con- 
venience, as 6c is not known a priori in 
AGN. Diffuse Galactic dust is consistent 
with a fraction al depletion of 6c = 0.6 
(IWhittetl 120031 ). Uncertainties about 6c 
will affect comparisons between columns 
derived from dust absorption and those 



^^Infrared rc-emission by meteoritic nanodiamonds 
is expect ed to generate emission band s at 3.43 and 



3.53 u m (|Van Kerckhoven et aLluOOa iJones et al 



inferred from the X-rays, as discussion in 
§ 15.31 and § 16.3.21 The dust-to-gas ratio 
arising from the C dust is 0.0031(5c (the 
solar neighborhood value assumi ng stan- 
dard ISM extinction is ~ 0.009, IWhittetl 



20031). 



20041 ) . Using the sp ectrum of 3C298 as a test 



spec 

case, Ide Diego et al.l (|2007i ) found no evidence of 



this emission. 



Matching the break-corrected UV 
segment with the X-ray SED. For 

each quasar, the column Nf}"- was var- 
ied, until the far-UV break could be repro- 
duced as closely as possible. Even though 
the curve D3 presents a steep decline red- 
ward of the cross-section peak (Fig. [6]), the 
extinction is not totally negligible in the 
near-UV. To preserve a good fit, one needs 
to slightly increase the index /5nuv with 
respect to the value of the dust-free case. 
Towards the optical domain, however, the 
extinction by nanodiamonds rapidly be- 
comes negligible. 

In Fig.[71 the thick black line represents 
the intrinsic (dereddened) SED, separately 
inferred for each quasar. The 'fluxed' X- 
ray models are the same as before and 
are represented by the green dot-dashed 
lines. To facilitate comparisons, the dust- 
free broken powerlaws of Fig.H] are re- 
peated (red dashed lines) . In TableH] we 
list the values of the intrinsic /3nuv index 
and other quantities related to the break- 
corrected SED. All the inferred quantities 
share the same superscript 'dia'. 

5.3. Discussion II: the break-corrected 
case 

After fitting the UV-break with dust- 
absorbed SEDs, we find that only 4 out of 
11 quasars are case A (a, e, f, g). There- 
fore, for a majority of objects, the adopted 
rollover function Cj-ow, with the same 6 = 
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—1.6 as in BOS, results in an extrapolated 
flux that is incompatible with the X-ray 
data. The possibility that the X-ray con- 
tinuum is significantly contaminated by an 
unresolved jet component would simply re- 
inforce this conclusion. Either the rollover 
must steepen further than assumed above, 
or a different functional form for the cut-off 
should be considered, such as an exponen- 
tial. For each object found to be caseB, 
we determined the steepening, (J^^^., that 
would revert it to case A. The values are 
listed in Col. 8 of TableH] and correspond 
only to upper limits, since further reduc- 
tions of 6 would still result in caseA's. 
The mean value of determined for the 
7 caseB's is —2.1, which is significantly 
steeper than assumed in BOS. 

In summary, the dominance of caseB 
objects in the break-corrected case indi- 
cates that the functional form adopted for 
the rollover should be revised. Indepen- 
dently of the adopted form, however, the 
flux recovery in the extreme-UV results 
in an ionizing photon luminosity that is 
higher than with the dust-free broken pow- 
erlaw discussed in previous §|11 a property 
further discussed in §[71 

It is interesting to compare the nanodi- 
amond dust columns N'^"' in Table[5] with 
the gas columns, , intrinsic to each 
quasar and inferred from the X-rays (§ 13. 2p . 
In the X-ray domain, it makes relatively 
little difference whether the metals lie in 
the dust or consist of free atoms, at the 
spectral resolution provided by Chandra. 
If the assumed dust screen covers both 
the X-ray and the UV source, the crys- 
talline dust column A?"^*" should be less or 
equal to the column of (intrinsic) gas in- 
ferred from the X-rays. A comparison of 



the two columns show that the condition 
A^^ > A^f " is satisfied for the 7 objects 
for which A^^ has been measured. It is 
also in agreement with the 2a upper lim- 
its of three other quasars (but not with 
the upper limit from quasar j). Both col- 
umn estimates rely on opacities computed 
for the solar metallicity case. Increasing 
or decreasing all the metals would not af- 
fect the column ratio Nf^°'/N-^. Reducing 
carbon depletion below the assumed value 
of 5c = 1, however, imphes proportionally 
larger Nf}°- columns. We find that a de- 
pletion in the range 0.1 < 5c < 0.3 allows 
the condition N^'^/N^ < 1 to be realized 
for the 7 objects for which A^^ has been 
measured. 

6. Combining the UV and X-ray 
segments assuming a universal 
SED 

6.1. The relevance of testing dust- 
reddened SEDs 

Photoionization requires a harder 
UV continuum. In their comprehensive 
review of accretion disks as the ultimate 
source of the AGN optic al-UV continuum, 
Koratkar fc BlaesI (119991 ) noted that 'bare' 
accretion disks, which are sufficiently hot 
to produce a significant fraction of hard 
ionizing photons, generally have an opti- 
cal/UV distribution that is too blue with 
respect to what is observed. These au- 
thors pointed out that the problem of disk 
models, however, does not lie in being able 
to fit the "red" optical/UV continuum, 
but rather in how to explain the extreme- 
UV emission, needed to power the high 



exci tation emission lines (see § 11.11 

and Korista et al. 1997bl ). To our knowl- 
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edge, there are no bare disk models that 
solve both problems simultaneously. This 
is the starting point for an exploration of 
whether residual dust extinction might be 
present in almost all luminous quasars, 
making their SED appear softer than it re- 
ally is. This would allow the fitting of bluer 
accretion disk models, which are more ef- 
ficient in generating a luminous and hard 
ionizing continuum. 



malization, since the grain composition 
and metallicity are different. The result- 
ing curve is shown as a continuous line in 
Fig. El For comparison purposes, we also 
show a silicate (MgFeSi04) grain model 
(long dashed-line) normalized to solar Si 
metallicit}0, but with otherwise the same 
grain size limits as the Pei model (i.e. 50 
- 2500 A). The two curves are similar in 
shape. 



SMC-type dust favored for SDSS 

AGN. An interesting result from the ex- 
tensive studies of AG N fro m the SDSS by 



Hopkins et al.l (120041 ) and iRichards et al. 



(120031 ) is that the optical properties of 
the dust in AGN, as inferred from op- 
tical and UV colors, are more akin to 
the properties found in the SMC than to 
the Ga l axMj^ . The recent simulations of 
WillottI (120051 ) supports the se findings. We 



adopt ed the conclusions of iHopkins et al. 



(120041 ) and calculated a SMC-type extinc- 
tion curve that is based on amorphous car- 
bon (AC) dust grains. The model consid- 
ers the grains to be spherical, with a size 
distribution oc a~^'^ within the range of 
50 < a < 1400 A. The complex refraction 
i ndice s n + ik are from iRouleau fc Martin 
(jl99ll ) for the AC type. We assumed so- 
lar metallicity for carbon and, to be defi- 
nite, full depletion of carbon onto dust (i.e. 
Sc = 1.0, see § 15.21) . The dust-to-gas ratio 
arising from AC dust is ~ 0.00345c- We 
fine-tuned the size range until our curve 
closely repro duced the shape of the SMC 
extinction of iPeil (119921 ). but not its nor- 



^^Thc extinction curve of Galactic dust is much shal- 
lower than that of the SMC and is characterized 
by a prominent 2175 A absorption bump, absent 
from the SMC curve. 



6.2. The hypothesis of a hard en- 
ergy distribution in the UV 

Assuming SMC-type extinction and a 
universal SED that is harder in the UV 
than what is directly observed, we proceed 
to investigate the amount of AC dust re- 
quired to provide an acceptable fit of each 
observed SED. We will keep the above 
nanodiamond dust component, in order to 
maintain a satisfactory fit to the UV break. 
More specifically, the working hypotheses 
behind the proposed test are the following: 

i- A pseudo universal index value for 
/3nuv- In their study of radio-quiet 



and radio- loud AGN, Gaskell et al. 



(120041 ) found what appears to be a 
low end cut-off in the distribution 
of near-UV indices. The hardest in- 
dices found corresponded to a^uv — 
—0.45. This value is quite similar 



"The dust model of |Pel (|l992l ) employed silicate 
grains but required 2-3 times higher Si abundance 
than available in the SMC interstellar medium. 
This one reason why we consider more attrac- 
tive the amorphous carbon grain composition. We 
chose to normalize our dust models to solar C 
abundances (C/H= 3.6 x lO"**) in order to be 
consistent with the metal content assume d for the 
opacities ( Morrison k, McCammoniri983 ) used in 
the X-ray model fitting f? [3?2)) . 
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to the optical-t o-UV index of the 



SPSS composite (jVanden Berk et al. 



2OOII). We define our Model I as a 



powerlaw with such an index (/^nuv ~ 
0.55), multiphed by the function 
Cfou, for consistency with previ- 
ous analysis. Five quasars with 
a measured index already harder 
than /Snuv — 0-^^ could not be in- 
cluded in this test and were left out. 
The SED from Model I is plotted in 
Fig. [8] (continuous line). The dot- 
ted line represents the same SED ab- 
sorbed by a crystalline dust screen 
with A^f'^ = 0.8 X 10^0 cm-2. We 
also defined a Model II based on 
an even harder SED. B05 encoun- 
tered a non-negligible fraction of 
quasars with a much hardei0 con- 
tin uum than the lirn it proposed 
by iGaskell et al.l (120041 ). some even 
harder than /3nuv = 0.8. Using the 



SPSS AGN sample, IVanden Berk et al, 
( I2OOII ) derived four composite SEDs 
of increasing hardness in the optical- 
UV. The hardest SED has /?nuv = 
0.75. We define Model II as a pow- 
erlaw with a similarly hard index 
of /^Nuv = 0.8, multiplied by func- 
tion Crou. It is plotted in Fig.O 
(continuous line) along with the 
popular AGN SEP publi shed by 
Mathews fc Ferlandl fll987f ). The 



dotted line represents Model II af- 
ter absorption by a crystalline dust 
screen of thickness A''^*" = 0.8 x 
1020 cm-2. 



-'^^BOS found a mean index for Class A quasars (see 
footnote[5]) of aNuv — —0.44 with a dispersion of 
0.21. 



ii- Determination of the AC dust col- 
umn. Initially, we adopt for each 
quasar the same column N'^"' of crys- 
talline dust required to reproduce the 
break (§[5]) and then determine the 
additional column of SMC-type dust 
needed for the fixed SED Model I 
or II (with either /5nuv = 0.55 
or 0.8, respectively) to fit each ob- 
served spectrum satisfactorily. For 
two objects (e and g) and only in 
the case of Model II, did we find that 
the fit to the far-UV break was af- 
fected when SMC-like extinction was 
included. For both, a small reduc- 
tion in crystalline dust column, iV]^*", 
restored the quality of the fit blue- 
ward of the break. For all the 9 other 
quasars, the column A'^^*'* was frozen 
to the values inferred from the break- 
corrected case of §|5l 

6.3. Discussion III: SMC-like ex- 
tinction in quasars 

6.3.1. Results from dust absorbed SEBs 

An interesting result is that, for both 
SEDs I or II, a combination of AC and crys- 
talline dust provides a fit to the UV spec- 
tra as satisfactory as that from a softer 
SED with nanodiamond dust only. This is 
illustrated by the three quasars a, f and 
g in Fig.[TOl The dashed lines represent 
amorphous carbon + crystalline dust ex- 
tinction, assuming the 'universal' Model II 
SED, while the dotted line represents ab- 
sorption by crystalline dust only, assuming 
the SED index /?nuv ^^^^ best fit each indi- 
vidual object (Col. 6 of TableH]). The thin 
long dashed line represents the intrinsic 
Model II SED (before extinction) for quasar 
c (3C351). Neither dust model (with and 
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without AC dust) results in a perfect fit of 
tlie spectra. However, taking into account 
tliat tlie true continuum level beneath the 
strong emission lines cannot be uniquely 
defined, we consider that both models, 
with and without AC dust, are equally suc- 
cessful. The kink around 1770 A in the AC 
extinction curve (Fig. [6]) results in a notice- 
able shoulder shortward of 1750 A in the 
transmitted flux of 3C454.3. The data are 
not inconsistent with such a feature. We 
conclude that accretion disk SEDs signifi- 
cantly bluer than observed are an enticing 
possibility. 

Since two SEDs of quite different hard- 
ness (I, II) can be fitted by simply vary- 
ing the amount of AC dust, we conclude 
that such a procedure does not lead to a 
unique description of the intrinsic SED. In 
Cols. 6 and 9 of TableO we list the to- 
tal^ absorption columns Njj and NjJ re- 
quired for a proper fit of each quasar spec- 
trum, assuming either the SED Model I or 
II, respectively. The corresponding color 
excesses E(B— V) due to intrinsic dust are 
shown in Cols. 7 and 10. We note that our 
E(B— V) values are relatively small, reach- 
ing at most 0.085 (object g), and that they 
are generally smaller than the Galactic val- 
ues (Col. 3). Nevertheless, such small color 
excesses strongly affect the f ar-UV SEP, 



as pre viously pointed out by iTripp et al. 
(I1994J ). who used an heuristic extinction 



curve that resembles the SMC extinction 
curve to determine upper limits on the 
reddening of 7 quasars. The Cols. 8 and 
11 list the column ratios of the absorb- 
ing material with respect to the one deter- 
mined in the X-ray domain. In Fig. [71 the 



thick cyan and yellow lines represent the 
intrinsic (dereddened) SED corresponding 
to Models I and II, respectively. 

SBG05 also explored the possibility of 
dust extinction by SMC- type dust. They 
found that the required extinction cannot 
remove the UV break without introduc- 
ing a large scale curvature to the dered- 
dened optical-to-UV region. They report 
that only 2 of their 17 AGN showed unam- 
biguous evidence of dust extinction. In our 
study, the break is specifically accounted 
for by crystalline dust absorption. When 
additional AC extinction is added, it does 
introduce some curvature, but not to the 
extent that it rules out AC extinction al- 
together. Presumably, the amount of AC 
dust required for our test is less than the 
amount probed by SBG05. 

As for the quasar sample studied by 
TZ02, the near-UV domain between 1300 
and 2000 A is characterized by a wide range 
of spectral indices among quasars. The dis- 
persion is as large as 0.57 about the mean 
value of /3nuv = 0.31. It would be inter- 
esting to explore whether variations in the 
extinction by amorphous carbon might not 
account for most of these variations in SED 
softness in the TZ02 sample. 

In their study of red and reddened 



uasars from the SDSS, Vanden Berk et al 



^^These columns correspond to the sum of the nan- 
odiamond and AC dust columns. 



2OOII ) derived one composite SED for sig- 
nificantly reddened AGN and four presum- 
ably unreddened composites that differ by 
the hardness of their continuum. These 
composites are characterized by optical- 
UV indices (p^) of 0.75, 0.59, 0.46 and 
0.24. The above authors did not rule out 
that dust extinction may have an impor- 
tant role in explaining the observed vari- 
ations in among their four composites. 
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However, they favor the idea that the dif- 
ferences are produced by changes in the 
intrinsic continuum, based on evidence of 
emission-hne trends with color. 

6.3.2. Requirements on the rollover shape 
implied by Models I and II 

We find that the extrapolation of the 
Model! SED results in case A in only 3 ob- 
jects (a, e, f) out of the 6 quasars for which 
we could carry out the test. In the case 
of Model II, however, all 11 quasars turned 
out to be case B. The conclusion therefore 
is that, while amorphous carbon extinc- 
tion favors the existence of an intrinsically 
harder SED, it also points to the need of 
a steep SED decline in the extreme-UV, 
much steeper than the one proposed by 
BOS. Alternatively, we may consider an 
exponential turnover as a more appropri- 
ate solution. Since no sharp bend is ob- 
served in the s pectra of the two qu asars 
HS 1700+6416 fiFechner et all l2006h and 
HE 2347-4342 fizheng et al.l \2004 \ down 
to ~ 250 A, such a sharp turnover can only 
take place shortward of this value. Since 
bare accretion disks predict a significant 
opaci ty above the ionizing threshold of He"*" 
(e.g. iHubeny et al.ll2000l ). we may reason- 
ably expect the proposed thermal turnover 
to occur above the He"*" edge (< 228 A). 

The ratio of the total dust column (crys- 
talline + amorphous) to the column in- 
ferred from the X-rays is given in Col. 8 
and 11 of Table[5l for Model I and II, re- 
spectively. In the case of Model II, this 
ratio is smaller than unity for the 7 ob- 
jects with measured absorption columns 
A^^. For the remaining 4 quasars, we 
have only upper limits for N§ . Note that 
each of them lies among the 5 objects with 



the smallest accumulated photon counts in 
the X-rays. Only for the object j does 
the lower limit column ratio exceed unity 
(by a factor of about two). The inferred 
Njj/N^ ratio can also be viewed as the 
minimum value of 5c required to equalize 
both columns. We conclude that, except 
for the 4 objects with upper limits, the 
absorption columns due to dust are over- 
all smaller than those inferred from the X- 
rays, provided the fractional depletion of 
C is not far from unity. This is consistent 
with but does not prove necessarily the 
concept that the same dust+gas compo- 
nent covers both continuum emission do- 
mains. Based on a s t udy o f BAL quasars, 
Gallagher fc EverettI (120071 ) have recently 
argued that the gas in the line-of-sight of 
the (inner) X-ray source is not the same 
than the one along the line-of-sight to the 
UV emitting region. 

Since the aox index is sensitive to ex- 
tinction of the UV-flux at 2500 A, we might 
expect that it correlates with the dust col- 
umn. In Fig.lTT], we plot the column A^^^ 
versus aox- There is no evidence of a 
correlation. Possibly, the intrinsic (dered- 
dened) aox varies significantly from object 
to object, which would tend to mask such 
a correlation. Assuming an aox intrinsic 
value of —1.6, the continuous line illus- 
trates the expected behavior of aox with 
increasing dust column. The slope is quite 
steep, which indicates that relatively mod- 
est variations in aox should be expected in 
any case. 

7. Comparison of photoionization 
models for the extended NLR 

In order to quantify the differences be- 
tween Models I and II and the TZ02 com- 
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posite, we computed photoionization cal- 
culations for the extended NLR (hereafter 
E-NLR). We favor the low-density regime 
(< lO^cm"'^), because the physics of the 
emission processes and radiation transfer 
are simpler and because it is unclear to us 
where to position the BELR with respect 
to the crystalline dust screen. The context 
and parameters of the models are described 
below. 

7.1. The E-NLR of radio-galaxies as 
a comparison reference 

The emission lines of the E-NLR, ob- 
served in projection on the sky or on the 
parent galaxy, are well suited for compar- 
ison with simple slab models. The rea- 
sons are twofold: i- in the case of the 
permitted lines, there is no need to per- 
form the difficult task of separating the 
NLR from the BELR component, ii- the 
gas densities of the emission regions seen 
in projection (outside the nucleus) are defi- 
nitely closer to the low-density regime than 
the nuclear NLR. The latt er is supported 



by the detailed analysis of iFerguson et al. 



( 119971 ) , who compared a wide range of line 
ratios from Seyfert I galaxies with those 
predicted by Locally Optimally emitting 
Cloud (LOG) models of the NLR. Their 
models are quite successful in reproduc- 
ing the high-excitation line strengths, the 
[O III] temperature and the correlation be- 
tween line widths and critical density ob- 
served in the (unresolved) NLR. As for 
the &NLR, the LOG scenario predicts 
[S ll] densities consistent with the low- 
density regime beyond ~ 200 pc and is 
able to reproduce the radial behavior of 
the [Oll]/[Olll] ratio in the E-NLR of 
NGG4151. These results indicate that spa- 



tially resolved emission-line components 
from the nuclear regions should be operat- 
ing in the low-density regime. 

The orient at ion- de pendent unifica tion 
scenario proposed by iBarthell (Il989l ). in 
which radio-galaxies and radio-loud quasars 
belong to the same parent population, 
justifies our proposed comparison of E- 
NLR calculations with the observations 
of the extended emission gas from radio- 
galaxies. This unification scenario has 
recently gained further support, follow- 
ing the mid and far-infrared c omparative 



studie s of F R2's and QSRs by iHaas et al. 



(I2OO4J . |2005| ). who used samples that en- 
compass a similar redshift and luminosity 
range for both classes. 

We assume photoionization by the nu- 
clear UV source as the main excitation 
mechanism of the E-NLR. A potential 
problem is that shock excitation might be 
contr ibuting significantly to the emission 
lines ( iDopita fc Sutherlandlll996l : lBinette et al, 



19851 ). How to determine the specific con- 



tribut ion fr o m sho cks is still an open ques- 
tion. iLaorl (119981 ) expressed strong reser- 
vations about the overall effi ci ency of fast 



shocks in AGN. iBest et al.l (l2000h 



pro- 
posed a scheme to separate the shock- 
excited case from the photoionized case, 
which is based on UV line ratios. In 
their proposed line-ratio diagram, the 
photoionized E-NLR of radio-galaxies oc- 
cupy a region characterized by the ra- 
tios Glll]/Gll]/ (1909A/2326A) > 3 and 
[Nelll]/[Nev] (3869A/3426A) < 2.4. In 
order to provide an observational dataset 
for the comparison of our calculations, 
we chose the spectrum from the z=2.36 
narrow-line radio-galaxy 4G— 00.54, which 
satisfies the first criterium (Glll]/Gll]/ = 
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8.1), so as to minimize the contribution 
from shocks. The hne ratios are hsted 
in Col. 3 of TableE] and the bibhography 
appears in footnote of the Table. We em- 
phasize that no attempt has been made to 
vary any of the model parameters in or- 
der to fit these ratios. Our main concern, 
above all, is to compare different SEDs with 
one another. 

7.2. Which SED to use: absorbed 
or unabsorbed? 

Since all quasars of redshifts > 1 found 
in the TZ02 sample show a far-UV break, 
and since the nanodiamond dust columns 
required to fit the break according to BOS 
cover a narrow rangj^. we infer that the 
covering factor of the dust is quite uni- 
form and near unity, that is, all line-of- 
sights are absorbed by similar amounts of 
dust. Hence, the SED to which the NLR 
of quasars is exposed must display a UV- 
break as well. For this reason, we will con- 
sider energy distributions from Models I 
and II that have already been absorbed 
by nanodiamonds (rather than the dered- 
dened form). Regarding the AC dust pos- 
tulated in §|6l since the observed UV spec- 
tral index varies considerably from quasar 
to quasar (TZ02), we infer that this hy- 
pothetical dust component is not uniform 
and presumably located further away from 
the ionizing source than the nanodiamond 
dust. In this case, it is reasonable to ex- 
pect that the NLR emission will be dom- 
inated by gas photoionized by radiation 
emerging from the line-of-sights that are 
least absorbed by the AC dust. For this 



reason, our E-NLR calculations will not 
consider SEDs absorbed by AC dust. The 
adopted energy distributions correspond- 
ing to Models I and II are plotted as dot- 
ted lines in Figs[8]and[9], respectively, while 
the SED that approximates the TZ02 com- 
posite is represented by the thick short- 
dashed line in Fig.O Both the TZ02 and 
the Model I SEDs include an X-ray com- 
ponent consisting of a powerlaw of index 
/3x = 0, assuming an aox of —1.5. 

7.3. Density stratified slabs 

The results of the model ing of the NLR 



in Seyfert I galaxies by iFerguson et al. 



( 119971 ) suggest that the partially ionized 
emission zones contain significant amounts 
of internal dust, while the fully ionized 
regions contain very little. When the ion- 
ization parametecj at the face of the slab, 
Uo, is high, radiation pressure becomes sig- 
nificant, and even small amounts of dust 
will lead to an internal density stratifica- 
tion. We therefore allowed our slab model 
to contain small amounts of amorphous 
carbon dust, at a level corresponding to 
2% only (i.e. 6c = 0.02), so as t o remain 



consiste nt with the iFerguson et al.l (119971) 
results. JPopita et all (l2002f ): iGroves et al.l 



( l2004a| ]bl) showed that for high Uo values, 
the ensuing density stratification result in 
an emission line spectrum that is relatively 
insensitive to the particular adopted value. 
In this case, Uo ceases to be a determi- 
nant parameter. This property of strat- 
ified models is particularly useful in the 
context where we are foremost concerned 



^"Thc mean dust column derived by B05 for the 
dominant Class A quasars was N20 = 1.05, with 
a dispersion of 0.29. 



•^^The ionization parameter at the face of the slab 
is defined as Uo = (pK/cn^, where (Pr/c is the 
density of ionizing photons impinging on the slab 
and the total H density at the face of the slab. 
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with the comparison of different SEDs. To 
be definite, we chose Uo = 0.2 and a slab 
face density of = 100 cm~^. Finally, 
we adopt the set of solar a bundances from 
Anders fc Grevessd (119891 ) 



7.4. Comparison of the three SEDs 

To compute line ra tios, we have used 



the code mappings ic ( Binette et al.lll985 



Ferruit et al.lll997l ) and assumed a slab ge- 
ometry illuminated on one side. For each 
ionizing SED that we considered, we calcu- 
lated the local equilibrium ionization state 
of the gas and integrated the ionization 
structure inward until less than 1% of the 
gas becomes ionized. We integrated the 
volume force exerted by radiation pres- 
sure within the sl ab, assuming hydros tatic 
equilibrium a s in iDopita et al.l (120021 ) and 
Binette et al.l (119971 ). Cols 5 and 6 of Ta- 
bleEllist the calculated line ratios with re- 
spect to H/S from SEDs corresponding to 
Model! and II, absorbed by nanodiamond 
dust, while Col. 4 represents calculations 
with the TZ02 composite. 

Tableini shows that the high excitation 
lines are of similar strengths in the case of 
Model I and of the TZ02 composite. This is 
due to the similar hardness of both SEDs, as 
reflected by the quantity (pn/^KeU, which 
is the ratio of ionizing photons above 1 ryd 
to those above 4 ryd (see values at bottom 
of table). The SED from Model II is signif- 
icantly harder, however, which is reflected 
in the calculated ratios. 

For an equal continuum flux near 1100 A 
(the onset of the break), the ionizing pho- 
ton flux, (fn, for Model I and II, is 23% and 
52% higher, respectively, than that from 
the TZ02 composite. These gains are quite 
modest, however. If, on the other hand. 



the observed softer index of the compos- 
ite were due to AC dust absorption (§I6]), 
the total gain in ipn would be much higher, 
reaching values of 62% and 270%, respec- 
tively (including previous factors). In- 
creasing the turnover energy beyond 54 eV, 
as suggested in § 16.3.21 would result in an 
intrinsically harder SED, which would in- 
crease the luminosity of the high excita- 
tion lines such as HellA1640, OviA1035 
and C IV A A 1549, and further contribute to 
resolve the softness p roblem reported by 
Korista et al.l (Il997bl ). By how much the 
turnover energy should increase would de- 
pend on the precise shape adopted for the 
turnover, a question that lies beyond the 
scope of the current study. 

8. Conclusions 

Using archived data from Chandra and 
HST-FOS, we have derived the UV to X- 
ray SEDs of 11 quasars. The aim of our 
analysis has been to constrain the behav- 
ior of the SED within the domain that is 
not directly observable, in the extreme UV. 
We explored the possibility that crystalline 
dust and amorphous carbon dust may ac- 
count for the surprising softness of quasar 
SEDs. More specifically, we have reached 
the following conclusions: 

1- When we fit the observed UV spec- 
tra with a broken powerlaw directly 
(dust-free case), we do not find a sin- 
gle object for which the soft X-rays is 
the prolongation of the far-UV pow- 
erlaw with the same index. In 9 out 
of 11 quasars, the extrapolated far- 
UV powerlaw lies at a flux level be- 
low or near the one actually observed 
in the soft X-rays. We find evidence 
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of soft X-ray excess in 5 objectcj. 

2- We simplify the nano diamond dust 
absorption model proposed by BOS 
to reproduce the far-UV break. In- 
stead of using two kinds of crystalline 
carbon grains (terrestrial cubic and 
meteoritic), we obtain a similar fit 
by using a dust model consisting of 
cubic diamonds only, but where the 
grain sizes cover a larger range (3- 
200 A). If the fractional depletion 6c 
within the gas is > 0.1, the dust 
columns required by nanodiamond 
(iwst models are not excessive, as they 
fall below the cold gas columns in- 
ferred from the X-ray data. 

3- Assuming that the far-UV break is 
due to nanodiamond dust absorp- 
tion, we explored the possibility of a 
shallow far-UV turnover at 18.5 eV 
to connect the far-UV with the X- 
rays, as proposed originally by BOS. 
We find that for most objects the 
steepening must be more pronounced 
than proposed by BOS, that is, an 
index change as steep as < —2. 
Given that two well studied quasars, 
HS 1700+6416 and HE 2347-4342, 
do not show any sharp break down to 
2S0 A, we also consider that the pro- 
posed rollover should take place at a 
significantly higher energy (> SO eV) 
than proposed by BOS. 

4- We find that adding SMC-type ex- 
tinction to the extinction from nan- 
odiamond dust leads to a bluer SED 



These are Pks 1127-14, B2 0827+24, 
Pks 1354+19, 3C454.3 and Pks 1136-13. 



and to a fit of the BBB of qual- 
ity comparable to the one achieved 
with nanodiamond dust alone. This 
opens the possibility that the intrin- 
sic quasar SED is much harder and lu- 
minous in the near and far-UV than 
previously considered. 
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A. Remarks on individual objects 



Additional information and special considerations about each quasar are described as 
follow: 

a.- PKS1127-14 . Th e X-ray observation of this quasar does not show pile-up. 



Siemiginowska et al.l (l2002l ) detected for the first time in the X-ray image of this object 
the presence of jet emission. For the main source they found a photon index rx=1.19 ±0.02 
in the spectral region between 2 and lOkeV, consistent with the best-fit value reported in 
this work within a level of 2a confidence. 

b.- PKS0405-123. (ISpinrad et aP Il993l : iBahcall et aP Il993bf ) found the presence of 
a Lya absorption line in the UV spectrum of this object. They associated this absorber 
with the rich cluster of galaxies where the source is located. After the Lyman limit 
V > 3.3 1 0^^ Hz ) there is an unusual excess with respect to a simple powerlaw model. 



Lee et al.l (Il995l ) found that Comptonization of the accretion disk cannot account for this 



excess. In our analysis of this source we created an additional mask between 717— 882 A to 
avoid the segment covered by this emission (see 13.11) . We note, however, that this does not 
produce a considerable change in the FUV index. 

The X-ray of this source shows jet emission (removed from our analysis). The source pre- 
sented a level of pile-up of 16.7%. We corrected the spectrum with the procedure described 
in §1221 

c- 3C 351 . This quasar shows the presence of Civ and Lya absorption lines (se e 



Fig.E] c) with 1700 k m s~ \ and 2200 km s ~^ bl ueshift, respectively f lBahcall et al.l Il993al ). 



Mathur et al.l (jl994al ) and iNicastro et al.l (119991 ) find, analyzing ROSAT data, the presence 



of an associated warm absorber in the X-ray band. This object prov ided the first suggestion 
for a common UV and X-ray absorber. Using a Chandra observation, iHardcastle et al.l (120021 ) 
fit the spectrum of this object with two power laws, attenuated by Galactic and intrinsic 
cold absorption. IHardcastle et al. (I2OO2I ) rejects a model consisting of ionized absorption 
that included only bound- free transitions (the model absori in XSPEC). One of the power 
laws in their model is associated with the nuclear quasar emission, while the other is assumed 
to be related to the jet emission. 

We analyzed the X-ray spectrum of this object using the approach described in § 12.21 The 
X-ray observation has ~ 11% of pile-up. Three X-ray prominent hot spots at 25 arcsec 
north-east from the main source were also detected. Those were excluded from our analysis 
to avoid contamination while getting the X-ray spectrum. 

The spectrum of this object, showed the presence of two data points with deviations larger 
than 2a near 2.5 keV. We did not included these points from our initial fit between 2.5-6 keV 
as they would change the measured value of /3x- Therefore, for this object, we fit only the 
range 2.51-6 keV. In this initial fit (that consists of a single powerlaw attenuated by galactic 
ab sorption) we found a spe ctral index /3x=0.6, consistent (within 2a) with the index reported 
by IHardcastle et al.l (120021 ). The extrapolation of this fit to the whole spectral range (0.3 
to 6 keV) showed that the spectrum is not well described with a single powerlaw. Strong 
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negative residuals between 0.5 <E < 2 keV are present, indicating the possible presence 
of ionized absorption, as previously suggested. The residuals further revealed the presence 
of an emission line consistent with the FeKa line. We thus fitted the spectrum with a 
powerlaw plus a gaussian (to account for the FeKa line) attenuated by both the WA and 
Galactic absorption. The WA was modeled using the code PHASE (IKrongold et al.ll2003l ) 
that includes absorption by both bound-bound and bound-free transitions. A single warm 
absorber could not account for the negative residuals, thus the warm absorber was modeled 
with two absorption compo nents. This is consistent with the findings for Seyfert galaxies 
(e.g. IKrongold et~aD l2nn5bl lal) . The best fit parameters for our fit are reported in TableO 
These model (d.o.f./x^ = 224/214) showed to be better than a model consisting in two power 
laws with intrins i c cold absorption (d.o.f./x^ = 227/255). 



Hardcastle et al.l (120021 ) reported the possible prese nce of an emission feature near 3.2 keV. 



We also notice the presence of this feature, and like iHardcastle et al.l (120021 ) could not con- 
strain its properties with a gaussian line. As our main interest is to constrain the con- 
tinuum properties, no further effort was made in explaining this feature [see discussion by 



Hardcastle et al.l (120021 )]. 



The above model for this object was carried out in the range between 0.32 to 6keV. The 
region between 0.24 and 0.32 keV showed several data points above of the power-law level 
that could not be fit by any model with a reasonable physical meaning. Furthermore, most of 
these points are located in the spectral region below 0.3 keV, where the spectral calibration 
for Chandra is not accurate. Therefore these points were excluded from our analysis. 

d. - 3C334. This quasar does not present pileup or jet emission in the X-ray domain. 

e. - B2 0827-|-24. The predefined mask in the UV domain (§ 13. ip was modified for this 
quasar. The spectral region between Lya and OviA1035, i.e., between 900-980 A was in- 
cluded to obtain a reliable FUV fit. It was also necessary to ignore the following regions: 
926.61-974.4, 896.5 897.5, 1090 to 1110, 1140.5 to 1142, 1310 to 1316 and 1440 to 1446 A due 
to the presence of absorption lines. Without these modifications to the mask, a harder FUV 
index (/9fuv= —1-01) would be obtained in the fits (we report an index of /5fuv= —1-27 in 
Table[2]). The X-ray observations of this object show the presence of a jet (removed from 
our analysis), and no pile- up. 

f. - Pks 1354-|-19. Also called 4C 19.44, this source did not presented pile- up in the X-ray 



spectr um. Jet emission was removed to avoid any source of contamination. iGambill et al. 
( 120031 ) analyzed the same observation and found variability by a decreasing factor of 3 during 
20 years (elapsed between the Chandra and the ROSAT observations). They reported an 
aox= —1.38 and photon indices of -1-1.39 and -1-1.88 for the hard and the soft spectrum, 
respectively. We found values consistent with those within a 2a level of confidence. 

g.- 3C454.3 An additional mask was introduced in the UV spectrum of this quasar 
between 877-885 A, to remove the absorption features present at a level exceeding 2a, and 
arising from the intergalactic medium. The X-ray spectrum presented 24% of pile-up, in 
spite the object was observed in use of the sub-array 1/4. The object further shows X-ray 
jet emission. 
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h. - OI 363 An additional mask was introduced in the range between 801-824 A of the UV 
spectrum to remove absorption features that were pr esent at a level exceeding 2a , and arising 
from the intergalactic medium. In the X-ray region, ISiemiginowska et al.l (120031 ) reported no 
pileup. However, with the pileup definition of this work, we found 16.3% of pileup on the 
same observation. The source presented evidence of jet emission. 

i. - Pks 1136— 13 The UV spectrum presented emission lines by Hell and Olll] in the 
range between 1640-1664 A. These were removed from the analysis in order to obtain a 
reliable fit of the continuum. In the X-ray band, the source does not show pile-up, but jet 
emission was found. 

PG 1 634-1 - 706 The X-ray observ ation does not present pile-up or jet emission. 



Page et al.l ( 12004J ): iPiconcelli et al.l (l2005l ) analyzed an XMM spectrum of this source be- 
tween 0.2 and lOkeV and found a photon index ~ 2 .19 consistent with the one reported by 
us within a level 2a of confidence. iPage et al.l (120041 ) further reported the presence of a soft 
excess, not required by the Chandra data analyzed here. 

k.- PG 1115-|-080 The UV spectrum of this source presents broad absorption lines. We 
then modify the UV mask to ignore these features in the spectral regions 886 to 892, 1097 
to 1105, 1138 to 114 and 1147 to 1149 A. In the case of considering those spectral ranges a 
flatter FUV index would be obtained (/3fuv~ —0.75, compared to the one reported in Table 
l2]/^Fuv= —0.65). In the X-ray region, the observation did not show pile-up or jet emission. 
W e note that in the X -ray band this quasar has two close sources [sources B and C in Fig. 4 
of lYoung et al.l (Il98ll )] due to gravitational lensing. These two sources were not present in 
the slit aperture of the HST-FOS observation. To be consistent with the UV data, sources 
B and C were removed from our analysis. 

This is the only quasar in our sample that required a broken powerlaw to fit the data. We 
first fit a single powerlaw between 2.5 to 6keV (rest-frame) absorbed by Galactic neutral 
gas, and obtain a spectral index /5x=0.86. However, the extrapolation of this model above 2 
keV resulted in a poor fit (d.o.f./x^ = 72/146). The best continuum model resulted from a 
fit with a broken powerlaw attenuated by Galactic absorption (d.o.f./x^ = 70/49). The best 
fit parameters for this model are: (1) an index /?x=0.86 in the 0.3 <E< 1.69 keV range (the 
spectral index was not changed from the previous value); (2) an index Pxhard = — 0-05^o;25 
in the 1.69 <E< 6keV range; and (3) an energy-break at E=1.69 keV. 
PI05 analyzed an XMM-Newton observation of this source. They used a single powerlaw 
and obtained a spectral index j3x= 0.15 (rx=1.85). The difference in the models arises from 
that we first performed the analysis in the soft X-ray region between 0.92 and 2.2 keV, and 
then extrapolated our model to the whole range of 0.3-6 keV (see § 13.21) . This is clear from 
middle panel of Fig. [3] A; where the residuals show a break around 2keV. 
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Fig. 1. — Comparison of different SED components in the UV and the X-rays: y-axis I'Fiy, 
lower axis: wavelength (A), top axis energy (eV). AH SEDs have been normalized to unity 
at Xnorm. = 2000 A. The SDSS AGN composite (gray solid line) was multiplied by 0.2 to 
avoid overcrowding. The red line represents the composite quasar SED from TZ02 that 
combines RQQ and RLQ quasars. The dashed purple line represents a broken powerlaw 
approximation of the TZ02 composite, using index values of 0.31 and —0.76, for 0mw and 
/ 3fuv, respectively. We also overlay the following theoretical SEDs: the iMathews fc Ferland 
fll987h SED (yellow dashed line). Models I and II absorbed by crystalline dust and introduced 
in §0 (dotted black lines), the thermal-like SEDs from K097 with an exponential cut-off at 
kT = 43 and 86 eV (short and long-dashed cyan lines, respectively). To avoid cluttering, 
UV SEDs end in a pointed arrow before the X-rays. The gray shaded area represents the 
wavelength region where scarce information exists. We can imagine different ways by which 
the extreme- UV connects within this data gap with the soft X-ray component. For a given 
X-ray SED, the index aox can be used to determine its flux level with respect to the UV 
domain. The thick orange markers at 2 keV delineate X-ray fluxes corresponding to arbitrary 
aox values of —1.3, —1.45 and —1.5, defined with respect to the 2500 A flux of the TZ02 
composite. As examples of AGN SEDs in the X-ray domain, we show two popular model 
fits: a broken-power law (navy blue line) and a blackbo dy-|-powerlaw (gre e n line ) fit. They 
correspond to unweighed averages of fits pe3&rmed by Brocksopp et al. ( 20061 ) and PI05, 
respectively (see § 11.21) . The blackbody+powerlaw X-ray model is favored in a larger number 
of AGN (PI05), and is adopted throughout this Paper. 




Fig. 2. — Rest-frame ultraviolet HST-FOS spectra of quasars a-k in uFjy as a function of 
frequency (bottom axis) or wavelength (top axis). Error bars are superimposed. The data 
have been rebinned to group 10 points per resolution element together. Separate powerlaw 
fits to the near-UV and far-UV regions are overlaid to the spectra (continuous lines). The 
two spectral indices /3^uv and /Sfuv a-iid the position of the UV break (A) are shown (see 
also Table E]). 
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Fig. 3. — Observer-frame X-ray spectra of the previous 11 quasars. The data are in 
counts/sec/keV. Top subpanel: final best-fit model superimposed to the observed count 
rate, middle subpanel: residuals from the initial model to the 2.5-6 keV region consisting in 
a powerlaw attenuated by Galactic absorption; bottom subpanel: residuals from the multi- 
component best-fit model. This Figure is available in color in the electronic version of this 
paper. 
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Fig. 4. — Overlay in uF^, vs. of the UV and X-ray spectral fits that were separately shown 
in Figj2]and FigjSl respectively. Black arrows at the top of each panel indicate the wavelength 
coverage of the original data. Quasar environment is assumed dust-free. Powerlaws of indices 
Pnvy^ Pf\jv (equivalent to actual slopes in uF^, plots) are used to describe the near- 

UV, far-UV and hard X-ray segments, respectively. A few representative wavelength and 
energy markers are shown above the top axisg^or both the UV and X-ray segments, dashed 
lines represent the best-fit models. Converging dotted lines indicate variations about these 
fits, assuming a 2a uncertainty for the input parameters, except for Xfj^, which was held 
constant. The symbol '@' denotes the UV break position, X'^l, and the circle the energy. 
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Fig. 5. — Rest-frame spectrum in Fx of PG 1148+549 (from BOS), multiplied by the scaling 
factor 0.45 x 10^"^ erg~^ cm^ s A. Pointers indicate the position of relevant emission lines above 
the continuum. The dotted line represents the intrinsic SED, which consists of a powerlaw 
with /3^{jy= 0.8, multiplied by the rollover function Cj-ou as defined in § 15.21 (Coincident ally, 
it is equivalent to Model II defined in § 16.21 ) Its flux is unity at the Lyman limit. The 
black continuous line represents the above SED absorbedhj nanodiamond dust, assuming the 
extinction curve D3 shown in Fig.[6]with N2o=^-0, while the dashed gray line represents a 
model with A^2o=0.8 that combines the Al and Dl curves (40% and 60%, respectively), as 
in B05. 
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Fig. 6. — Cross-sections (normalized to H) of different dust models as a function of wave- 
length. Extinction curves labeled Al (meteoritic nanodiamond) and Dl (terrestrial cubic 
diamond) were used by BOS to reproduce the far-UV break of quasars. Curve D3 (continu- 
ous black line) is used throughout this Paper to fit the quasar break near iiooA (see § 15.21 
and Fig.E]). It consists of terrestrial nanodiamonds that cover a wider range of 3-200 A as 
compared to the Dl curve (3-25 A). The long-dashe d lin e is an extinction model of the SMC, 
using silicate grains with the same size range as in iPeil (119921 ) . The normalization assumes 
solar metallicity of Si, and the resulting extinction curve has been scaled by a factor 2 for 
comparison purposes. We will adopt a different dust model that behaves similarly to that of 
the Pei model, but consists of amorphous carbon grains instead (solid line labeled AC). The 
normalization for AC assumes solar abundance of C and full depletion onto grains {6c = 1). 
The curve AC will be used to deredden the near- UV continua of quasars, as described in 
§E21 
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Fig. 7. — Overlay in vFj^ vs. z/ of the X-ray spectral segment (green dash-dotted line) with 
different UV models. The thick black line: overlay of the intrinsic UV SED of index /3nuv 
needed to fit the break assuming nanodiamond dust only (§[5]); the cyan line: overlay of the 
intrinsic UV SED Model! (/?nuv — 0.55) that fits the UV spectra when absorbed by both 
amorphous carbon and nanodiamond dust (§ 16.21) : the yellow line: overlay of the intrinsic 
UV SED Model II (/^Nuv ~ O-^) ^^^^ ^'^^ ^¥ spectra when absorbed by both amorphous 
carbon and nanodiamond dust (§ 16.21) . All UV SEDs include a rollover at 670 A defined by 
function Crou, which steepens the powerlaw index by amount —1.6 (§ 15.21) . Other symbols 
have the same meaning as in Fig. HI For comparison purposes, we also overlay the dust-free 
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Fig. 8. — Comparison of Model I with other spectral energy distributions. The continuous 
line represents the intrinsic SED from Modell, as defined in § 16. 2[ It consists of a powerlaw 
with /3nuv — 0.55, multiplied by the turnover function Crou defined in § 15. 2[ This turnover 
at 18.5 eV (4.48 x 10^^ Hz) consists of an index steepening (by 6 = —1.6), bringing the far- 
UV index progressively towards /3puv = —1-05. The dotted line represents the same SED, 
absorbed by nanodiamond dust only, assuming iV2o=0.8 and the extinction curve D3 (Fig.[6]). 
The long dashed line represents the softer of the two SEDs used by K097 in their BELR grid 
of models. It peaks at 22 eV. The short dashed line is a broken powerlaw approximation of 
the composite SED of TZ02 (with /?nuv = 0.31 and Pfuy = —0.76). In order to include the 
X-rays in the photoionization calculations of §[3, a truncated powerlaw of index f3^ = 0.0 
(Fx = —2) was appended to Modell and the TZ02 composite, so as to obtain an aox of 
— 1.5 in both cases. All distributions are normalized to unity at 2500 A. 
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Fig. 9. — Comparison of Model II with different spectral energy distributions. The contin- 
uous line represents the SED Model II, as defined in § 16. 2[ It consists of a powerlaw with 
/^Nuv ~ 0.80, multiplied by the same turnover function Crou as defined in § 15.21 The dot- 
ted line represents the same SED absorbed by nanodiamond dust only (extinction curve 
D3), assuming N-)n=0-8. The dot-dashed line line represents the AGN SED proposed by 



Mathews fc Ferlandl ( 19871 ). All distributions are normalized to unity at 2500 A. 
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Fig. 10. — Spectral energy distributions in Fx as a function of A for the quasars c: 3C351, /: 
Pks 1354+19 and g: 3C454. The dotted hues correspond to absorption by crystalhne carbon 
grains only, assuming an intrinsic powerlaw index /3nuv 0.28, 0.02 and —0.09, respectively. 
The thick short-dashed lines represent absorption models that incorporate amorphous as well 
as crystalline carbon grains. The SED corresponds to Model II for all three objects, that is a 
powerlaw of index /3^uv = 0.8 in the UV with a rollover at 670 A. The thin long-dashed line 
represents such an intrinsic SED, normalized to fit the spectrum of 3C351 when absorbed by 
dust. 



41 




I I I I I I I I I I I I I I I I I I I 

-2 -1.8 -1.6 -1.4 -1.2 -1 

«ox 

Fig. 11. — For each quasar, the absorption column inferred by fitting Model II, N^j ^ versus 
measured ctox- Using the absorption cross- section of AC, the continuous line gives the 
expected behavior of ctox if all objects had the same intrinsic value of —1.6. The dashed 
line illustrates the hypothetical case of a threefold increase of the dust column, assuming an 
intrinsic index of —1.8. 
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Table 1 

Observation Log for the 11 quasars 



Common 


Fig 


z 




rtaaio- 


Spectral 


Chandra 


Exp. Time 


Total 


Frame 


Chandra 


xlD 1 


Name 


ijabci 




lU cm 


loudness 


Ciass^ 


lU 


(ksec) 


Counts^ 


i ime 


Date 


Date 


(1) 






(.4) 


(.5) 


(,f>) 


(,') 




(.9) 


(10) 


/"1 1 ^ 
(.11) 


/"I OA 

(lij 


Pks 1127-14 


a 


1.18 


4.09 


RLQ 


A 


866 


27.3 


16370 


0.4 


2000-05-28 


1993-01-01 


Pks 0405-123 


b 


0.57 


3.8 


RLQ 


A 


2131 


8.6 


11994 


0.4 


2001-07-22 


1991-07-01 


3C351 


c 


0.3712 


2.49 


RLQ 


A 


2128 


50.92 


9450 


1.7 


2001-08-24 


1991-10-22 


3C334 


d 


0.56 


4.28 


RLQ 


A 


2097 


32.4 


8928 


0.85 


2001-08-22 


1991-09-07 


B2 0827+24 


e 


0.939 


3.62 


RLQ 


A 


3047 


18.2 


6914 


0.4 


2002-05-07 


1997-10-28 


Pks 1354+19 


f 


0.720 


2.23 


RLQ 


A 


2140 


9. 


6445 


0.4 


2001-01-08 


1992-02-26 


3C454.3 


9 


0.850 


6.41 


RLQ 


B 


3127 


4.9 


4716 


0.8 


2002-11-06 


1991-09-11 


OI 363 


h 


0.63 


4.18 


RLQ 


B 


377 


27.6 


4291 


3.2 


2000-10-10 


1997-05-15 


Pks 1136-13 


i 


0.560 


3.6 


RLQ 


A 


2138 


8.9 


3720 


0.4 


2000-11-30 


1992-01-30 


PG 1634+706 


j 


1.330 


4.54 


RQQ 


A 


1269 


10.8 


3113 


0.4 


1999-08-21 


1991-11-03 


PG 1115+080 


k 


1.718 


4.01 


RQQ 


A 


363 


26.4 


1814 


3.2 


2000-06-02 


1997-01-22 



''Webtool COLDEN |http:// cxc . harvard .edu/toolkit/ colden . j sp [ which is based on radio Himaps of lDickev fc LockmanI lll990l ). 
OJ ''Spectral class based on the far-UV break, following B05. 

'^Extracted from Sherpa's command "show" using a 3 arcsec circle region around the source. 

■^See table with CDD frame-time (sec) for standard Subarrays | http:/ /cxc. harvard. edu/proposer/POG/| 



Table 2 

Broken powerlaw fit to the observed UV continuum 



Common 


Object 


S/N^' 


Spectral 


\dfe 
^brk 


Near-UV'' 


FaJT-UV'' 


Ab 
A" 


Name 


Label 






Range 


(A) 






(10-15) 


(1) 


(2) 


(3) 




(4) 


(5) 


(6) 


(7) 


(8) 


Pks 1127-14 


a 


8, 


,5 


681-1498 


1022 


^U.iy+o.or, 


-1.82 + 0.08 


3.60 


Pks0405-123 


b 


18, 


,5 


733-2083 


1300 


0.74 + 0.01 


-0.40 + 0.03 


38.1 


3C351 


c 


10, 


,3 


840-2386 


1039 


0.28 + 0.01 


-1.26 + 0.20 


20.1 


3C334 


d 


8, 


,9 


884-3075 


1285 


0.29 + 0.01 


1 -0,24 
-•-+0.23 


6.55 


B2 0827+24 


e 


9, 


,5 


861-1690 


988 


-0.36 + 0.07 


-1.27 + 0.41 


4.87 


Pks 1354+19 


f 


9, 


,3 


788-2780 


1303 


0.02 + 0.02 


—0 70-0.16 
"■'"-HO. 14 


6.03 


3C454.3 


g 


11, 


,9 


877-2575 


1274 


-0.09 + 0.01 


1 79-0.27 
■^■'^-HO.25 


5.32 


OI 363 


h 


13, 


,6 


789-1406 


1097 


0.04 + 0.25 


-1.89 + 0.28 


8.73 


Pks 1136-13 


i 


12, 


,2 


815-2100 


1300 


0.78 + 0.02 


-0.95 + 0.14 


6.94 


PG 1634+706 


j 


12, 


,5 


809-1412 


1143 


0.40 + 0.13 


-1.36 + 0.11 


21.0 


PG 1115+080 


k 


15. 


,8 


817-1759 


1069 


0.34 + 0.01 


-0.65 + 0.02 


4.62 



■^Signal to Noise (F,. /error-bar) at 1300 A. 

''The powerlaw is normalized as follows: = A {v/Vj.^f)'^"''' = A {i'/i'ref)^^''~^\ where 
Uj-gf = c/1197 X 10-* cm. We recall that /3i, is defined in the i/Fi, plane, with /Si, = av + 1. 



Table 3 

X-RAY SPECTRAL FITS (OBSERVER-FRAME) 



Common 
Name 



Fig. 
Lab 



Best Fit 
Model Descr. 



PL Neu. Gas d.o.f/ x^'"^^ 



Ionized Gas 
LogU / N§/ z 



zBB 
kT/N(«) 



G 



fwhm/pos/A(-'') 



d.o.f./xVQ 



(1) 



(■■-!) 



{■>) 



(T) 



(11) 



Pks 1127-14 
Pks 0405-123 
3C351 

3C334 

B2 0827+24 

Pks 1354+19 

3C454.3 

OI363 

Pks 1136-13 

PG 1634+706 

PG 1115+080 



a GA*IA*(PL+zBB) 

b GA*IA*PL 

c GA*IA1*IA2*(PL+G) 

d GA*IA*PL 

c GA*IA*(PL+zBB) 

f GA*IA*(PL+zBB) 

g GA*(PL+zBB+G) 

h GA*IA*(PL) 

i GA*{PL+zBB+G) 

j GA*(PL+G) 

k GA*BP 



1.11 + 0.12 

n OO-0.66 

0.60 + 0.16 



0.35 ± 0.23 

0.71 + 0.22 
0.82 + 0.23 
1.32 + 0.40 
0.44 + 0.72 
0.93 + 0.36 
0.08T 



-0.39 
^+0.40 

0.86 ± 0.36 8 



4.02 
1.54 
2.8 

2.60 
2.95 
4.95 
1.89 
0.66 
2.30 
2.6 
1.67 



27.347 



-7.13 
*+8.8 

4.U9_,_4 23 



13.057 



'+1.53 
-5.52 
'+7.26 
fi 44-12.66 
"■^^+7.71 

< 43.38 

'•^^+6.05 

< 2.40 

< 1.16 

< 14.28 



285/608 
49/33 
232/900 

177/229 
172/219 
166/206 
67/117 
61/54 
220/396 
93/92 
71/77 



1.38 + 0.03/363.1 + 1.07/0.36 
0.90 ± 0.21/26.3 + 1.32/0.32 



0.63/2.44 



282/272/0.65 
48/28/0.99 
224/214/0.67 

176/149/0.93 

169/128/0.99 
163/143/0.86 
62/44/0.96 
60/46/0.76 

0.25/1.35 0.01/ 3.35/2.57 115/90/0.71 
••■ 0.001/2.84/28.92 90/88/0.54 
70/49/0.97 



0.46/1.13 
0.28/1.54 
0.43/1.49 



0.08/2.66/0.23 



0.42/2.14/1.20 



^Index of powcrlaw fits carried out, considering the 2.5-6 kcV (rest-frame) region only. We recall that /3 indices are defined in vFi, vs u, with /3x = ckx + 1 = 2 — Fx- 
''Powerlaw normalization defined at 1 keV in 10"* photons keV~^ cm,~^s~^ . 

^Column density of Neutral Gas x 10^" cm~^ due to intrinsic absorption and associated 2cr error. 

'^Statistics during initial fit; Power Law attenuated by Gala<;tic absorption and extrapolated between 0.3 to 6 keV (Observed Frame). 

■^Black Body Temperature (keV) and normalization in units of 10~^ Lgg/ D'^^, where L39 is the source luminosity in units of 10^^ ergs and Dio is the distance to the 
source in units of 10 kpc. 

'l-D unnormalized Gaussian function. Listed: full-width at half-maximum (keV), mean position (keV) and amplitude (10^*). 
^Broken Power Law. Spectral Index after breaJf (1.69 keV) is /3x= ~0-05+o 25 ^"'^ energy of reference at 0.5 keV. 



Table 4 

Parameters for the four different UV sed Models 



Quasars 




Dust-free Model 


Break-corrected Model 


AC Dust Modell 


\^c\Ae^ TT 


Common 


Object 




^NUV 


join 


adia 
^NUV 


join 


^dia 


JOtlX 


Tfi,0,X 


s:J 1 



Name 


Label 






eV 




eV 






eV 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


Pks 1127-14 


a 


-1.11 


-0.19 


146 


0.4 


375 


-1.21 


496 


-1..39 


-1.78 


Pks 0405-123 


b 


-1.78 


0.74 




0.8 




-2.73 


(a) 


(a) 


-2.86 


3C351 


c 


-1.56 


0.28 


432 


0.34 




-1.87 




-2.22 


-2.60 


3C334 


d 


-1.46 


0.29 


327 


0.6 




-1.85 


(a) 


(a) 


-2.15 


B2 0827+24 


e 


-1.28 


-0.36 


223 


0.1 


253 


-0.97 


583 


-1.60 


-2.00 


Pks 1354+19 


f 


-1.29 


0.02 


429 


0.2 


244 


-1.0 


482 


-1.54 


-1.90 


3C454.3 


g 


-1.31 


-0.09 


216 


0.1 


353 


-1.10 




-1.81 


-2.20 


OI363 


h 


-1.72 


0.04 


194 


0.29 




-1.92 




-2.32 


-2.72 


Pks 1136-13 


i 


-1.34 


0.78 


341 


0.78 




-1.83 


(a) 


(a) 


-1.89 


PG 1634+706 


j 


-1.64 


0.40 


247 


0.79 




-2.17 


(a) 


(a) 


-2.16 


PG 1115+080 


k 


-1.62 


0.34 




0.64 




-2.07 


(a) 


(a) 


-2.23 



''Model could not be applied to these quasars, since their near-UV index is already harder than the target value of /S^uv ~ 0.55. 



Table 5 

Absorption columns due to dust and metals 



Quasars 


Galactic Abs. 


X-rays 


ND Dust 


AC Dust Model I 




AC Dust Modem 


Object 


H 






n. 


All 
H 


tj[D-V ) 






tLyti- V ) 


\rll 1 atX 


Label 


1020 cm- 2 




1020 cm-2 


1020 cm-2 


lO^o cm-2 




1020 cm- 


-2 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


a 


4.09 


0.069 


27.3 


1.4 


2.6 


0.017 


0.095 


4.6 


0.045 


0.17 


b 


3.8 


0.064 


4.1 


1.0 








1.25 


3.6x10-3 


0.31 


c 


2.49 


0.042 


389 


0.60 


2.3 


0.024 


0.0059 


3.8 


0.045 


0.010 


d 


4.28 


0.073 


4.7 


1.5 








3.0 


0.021 


0.64 


e 


3.62 


0.061 


13.1 


0.90 


3.8 


0.041 


0.29 


5.4 


0.065 


0.41 


f 


2.23 


0.038 


6.4 


0.60 


3.5 


0.041 


0.54 


4.5 


0.055 


0.69 


g 


6.41 


0.11 


< 43.4 


0.70 


4.8 


0.053 


>0.11 


6.7 


0.085 


>0.15 


h 


4.18 


0.071 


7.2 


1.5 


3.5 


0.029 


0.48 


5.2 


0.054 


0.72 


i 


3.6 


0.061 


< 2.4 


1.2 








1.4 


2.3x10-3 


>0.58 


j 


4.54 


0.077 


< 1.2 


1.8 








2.1 


4.3x10-3 


>1.8 


k 


4.01 


0.068 


< 38.5 


0.90 








1.2 


4.3x10-3 


>0.031 



Table 6: Line ratios from photoionization calculations'^ 



Line ID 


A 


E-NLR^ 


TZ02 


Model P 


Model ir 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


LvfT 


1216 


7.1"^ 


12.3 


13.7 


13.8 


B.f3 


4861 


1.0 


1.00 


1.00 


1.00 


Ha 


6563 


2.9 


2.93 


2.92 


2.91 


Hen 


1640 


0.58 


1.2 


1.3 


1.7 


Hen 


4686 




0.15 


0.17 


0.22 


Cii] 


2326 


0.10 


0.05 


0.03 


0.05 


Cm] 


977 




0.029 


0.020 


0.042 


Cm] 


1909 


0.26 


0.5 


0.5 


0.8 


Civ 


1549 


0.79 


1.4 


1.3 


2.1 


[Niil 


6583 


0.83 


0.21 


0.18 


0.23 


Nlll] 


991 




0.0018 


0.0012 


0026 


Nlll] 


1749 




0.024 


0.023 


0.038 


Niv] 


1485 




0.12 


0.11 


0.18 


N V 


1240 


0.34 


0.36 


0.33 


0.54 


[0 il 


6300 


0.29 


0.19 


0.09 


0.16 


L J 


3727 


1.3 


0.56 


0.54 


0.69 


[Om] 


4363 




0.057 


0.060 


0.10 


[Oiii] 


5007 


9.8 


10 


11 


14 


Om] 


1665 


0.10 


0.084 


0.075 


0.15 


Ovi 


1035 




1.3 


1.1 


2.6 


[Sii] 


6723 


0.60 


0.36 


0.18 


0.27 


Mgll 


2800 




0.29 


0.22 


0.32 


Sim] 


1887 




0.0023 


0.0022 


0.0051 


[Ne III] 


3869 




0.57 


0.66 


0.87 


[Neiv] 


2424 


0.61 


0.2 


0.2 


0.27 


[Nev] 


3426 




0.61 


0.63 


0.93 


[Fevii] 


6087 




0.00085 


0.00089 


0.0011 


Fex] 


6374 




0.0076 


0.0046 


0.025 


V^Hell/V^H 






0.088 


0.093 


0.127 


(Te) 


(K) 




9450 


9300 


10100 


(He) 


(cm ^) 




595 


565 


560 
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Line ratios are expressed relative to H/3. All models are ionization bounded and assume 
a 'front' density = 10^ cm~^, solar metallicities and dust grains made of amorphous 
carbon corresponding to a depletion of only 6c = 0.02 (see § 17.21) . 

^ Line fluxes, relative to H /? , for the z = 2. 3 6 nar row- line radio-galaxy 4C— 00.54, taken 
from lVernet et al.l (120011 ): llwamuro et al.l (120031 ) and Humphrey et al. (in prep.). The 



relatively high-ionization state in the E-ELR of this source (Clll]/Cll] = 8.1) (see 
§ l7.ip . taken together with the unperturbed kinematics (FWHM~ 700kms~^ and large 
ra dio size (~ 200 kpc) , suggests that it is not strongly affected by shock-ionization 



(IHumphrey et al.ll2006l ). Since B.j3 was not detected in the spectrum of this source, we 
have assumed that Ha/ H/j = 2.9. The line ratios are corrected for Galactic reddening 
using E(B-V) = 0.026 (IVernet et al.ll200lh . 



Both SEDs labeled Model I and II have been absorbed by nanodiamond dust and are 
plotted as dotted lines in Figs.[H]and[9], respectively. 

The Lya is most likely absorbed, either by internal dust or possibly by an extended 
but local H I absorb er, which is a common phenomenon in high redshift radio-galaxies 



(IBinette et al.ll2006l . and references therein). 
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